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ABS1BACT 


This  thesis  reports  on  methods  for  solving  traveling 
salesman  problems  with  time-window  constraints-  Two  types  of 
time  windows  are  considered:  hard  time  windows,  which  are 
inviolable,  and  soft  time  windows,  which  are  violable  at  a 
cost.  For  both  cases,  we  develop  several  heuristic  proce- 
dures, including  some  that  are  based  on  Stewards  [Ref-6] 
effective  heuristics  for  the  traveling  salesman  problem 
without  time-window  constraints.  In  addition,  we  develop 
exact  algorithms  for  each  case,  which  are  based  on  the 
state-space  relaxation  dynamic  programming  metnod  of 
Christo fides,  Mingozzi,  and  Toth  [Ref.5].  Computational 
experience  is  reported  for  all  the  heuristics  and  algorithms 
we   develop. 
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I-     INTfiODDCTION 

A-       OVEBVIEW 

Consider  a  traveling  salesman  having  to  visit  n  cities 
or  customers.  He  starts  from  a  depot  and  needs  to  visit  each 
of  the  other  n-1  cities  only  once  and  then  return  to  the 
depot.  The  ccst  of  traveling  between  any  pair  of  cities 
(expressed  in  terms  of  distance,  time  or  cost,  etc),  say 
from  city  i  tc  j,  is  given  as  c?;  in  a  cost  matrix  C.  The 
problem  is  to  design  a  tour  through  the  n  cities  that  would 
minimize  the  total  cost  of  the  tour.  This  is  known  as  the 
Traveling  Salesman  Problem  which  is  a  well-known  classical 
operations    research    problem. 

The  TSP  is  called  Euclidean  when  the  cities  that  must  be 
visited  all  lie  on  the  same  plane  and  the  cost  of  traveling 
between  any  pair  of  cities  is  the  Euclidean  distance  between 
them. 

The  TSP  is  an  NP-complete  problem  [ Ref .  1r  2].  All  known 
exact  solution  methods  have  a  rate  of  growth  of  the  computa- 
tion time  which  is  exponential  in  n.  On  the  other  hand, 
heuristic  solution  methods  have  a  rate  of  growth  of  the 
computation  time  which  is  a  lew  order  polynomial  in  n  and 
have  been  experimentally  observed  to  perform  well.  For  this 
reason,  there  has  been  an  extensive  amount  of  research 
directed  at    TSP    heuristics. 

In  this  thesis  we  consider  adding  time  window 
constraints  to  the  TSP.  That  is,  if  t;  is  the  time  tnat  the 
salesman  visits  city  i,  then  t^must  satisfy  l'L  <  tL<  u^ , 
where  1-L  and  u«  are  the  specified  lower  and  upper  bounds  of 
a  time  window.  This  problem  is  not  as  well  studied  as  the 
unconstrainted  TSP,  but  there  have  teen  a  few  approaches 
used   on    the    problem. 


Psaraftis  [Ref.  3]  has  presented  a  dynamic  programming 
model  and  solution  procedure  for  two  dial-a-ride  problems, 
which  are  similar  to  time-window  constrained  TSPs.  Baker 
[Ref.  4]  has  presented  an  exact  algorithm  using  a  branch  and 
bound  procedure  which  is  effective  for  very  small  n. 
Christofides  et  al.  [Ref.  5]  have  presented  a  dynamic 
programming  state  space  relaxation  procedure  to  compute 
bounding  information    within    a   branch    and   bound   algorithm. 

The  objective  of  this  study  is  to  develop  exact  and 
heuristic  algorithms  which  will  provide  an  optimal  or  near 
optimal  tour  that  visits  each  city  in  its  given  time 
window.  He  are  given  a  depot  location,  a  set  of  x,y 
co-ordinates    for   n   cities   an d   a   set    of    time    windows   - 

A  common  application  of  the  TSP  is  in  veiiicle  routing 
problems.  A  set  of  customer  orders  must  be  partitioned 
among  several  vehicles.  Given  a  partition,  the  problem  tnen 
decomposes  into  one  TSP  for  each  vehicle.  Because  of  this 
prospective  application  and  in  deference  to  the  difficulty 
of  large  TSPs  (with  or  without  time  constraints),  we  confine 
our  research  and  computation  to  small-scaled  problems  (less 
than  3  0   nodes). 

We  consider  two  different  kinds  of  time  window 
ccnstraints;  hard  time  windows  and  soft  time  windows.  Hard 
windows  cannot  be  violated.  Soft  windows  can  ce  violated, 
but  a  penalty  cost  must  be  paid  when  they  are-  The  penalties 
can  be  defined  individually  for  each  customer,  and  they  can 
differ  for  early  and  late  arrivals.  Generally,  the  penalty 
for  arriving  before  the  lower  time  window  bound  is  much  less 
than  tne  penalty  for  arriving  after  the  upper  bound-  In 
Chapter  III,  we  present  the  hard  time  window  approach  and  in 
Chapter    IV,    we    present   the   soft   time    window    approach. 

We  developed  several  Fortran  programs  for  solving  the 
TSP  and  time-constrained  TSP.  For  the  TSP,  we  use  Stewart's 
[Ref.       b]   recent      heuristics,      CCCO    and   CCAO.         For    the    time 
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constrained  ISP  problems,  we  develop  some  new  heuristics, 
seme  of  which  are  modification  of  Stewart's  heuristics  for 
the  unconstrained  problem.  We  also  developed  exact  algo- 
rithms for  both  hard  and  soft  windows  using  Christofides  et 
al. 's  [Ref.  5]  method  of  state-space-relaxation  dynamic 
programming  and  branch  and  bound.  This  is  described  in 
Chapters  III    and  IV. 

Finally,  we  describe  a  hybrid  of  the  heuristic  and  exact 
programs.  Tie  hybrid  uses  the  overall  structure  of  the 
exact  program,  but  the  upper  bounds  are  obtained  with  the 
heuristic.       This   is   discussed    in   Chapters   III    and    IV. 

B.       THE    TfiAVEIING    SALESMAS    PBOBLEH 

A   tour    is      a   chain   which    passes    through   all      the   n    nodes 

and    in    wnich    the  first  and  the  last    nodes   coincide.         A    tour 

is   also   known   as   a    Hamiltonian  cycle. 

Let   a   tour   be   denoted    by    t  =    (i    ,    i    ,. .  . ,    i    ,    i    )       and 

12  n  1 

the    cost  of    this   tour    be 

n-1 
C  (t)     =        *•>""    c  +      c  . 

**-—>      i      i  i    i 

j  =  1         j      j+1  n    1 

Here         {   i    ,    i    ,...,i    )    is      a      permutation    of      the      inteaers 
12  n 

from    1     to  n,      giving      the      order      in      whica      the      cities   are 

visited. 

The  Traveling  Salesman  Problem  can  be  defined  as 
follows.  Given  a  graph  G  =  {  N,A  }  composed  of  a  set  of 
nodes  N,  a  set  of  arcs  A  connecting  these  nodes,  and  a  cost 
(distance  )  c>»  associated  with  each  arc  (i, j)  in  A.  The 
TSP  is  the  protlem  of  finding  the  minimum  cost  tour  of  the 
nodes  in  N.  The  following  mathematical  formulation  of  the 
TSP  is  from  Stewart  £Ref.  6]. 
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MIN  y>  C  X 

— !  ij      ij 


(1.1) 


1 

#  D 

S.T 

13 

n 

y-y..     - 

i>  - 1 

5i 

y . 
13 

-u    x           <      0 

j    =    1, .. . rn 

(1-2) 

i   =    1,...,n 

(1.3) 

i  =    1,.. .,n  (1.4) 


i   =    2, .- . , n  (1-5) 


j    =    1,-.. 

,,n 

i*j 

x      =   Or  1 
ij 

for    all      (i,  j) 

y      >   0 
ij 

for    all       (i,  j) 

(1-6) 
(1,7) 


where 


x        =  /  1  if    arc  (i,  j)    is   on    the    toar 

iJ 

\  0  otherwise 

y  are   continuous   variables    that    force   the    final 

ij 

solution    to   be  on    the      tour 

(   i.e.    include  every    node    in  the    same   route    ) 

u  is   a   number   >   n-  1    ,    and 

n  is   the    number    of    nodes    in    the  set    N. 

The  constraints  (1.2)  and  (1.3)  ensure  that  each  node 
will  be  visited  exactly  once,  while  constraints  (1.4)  and 
(1.5)  force  the  final  solution  to  be  a  single  tour  that 
starts  and  ends  at  node  1  (depot) .  This  formulation  is  not 
directly  used  in  our  TSP  programs,  but  is  of  interest  in  a 
general    discussion   of    the    problem. 
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C.       TSP    WITH    1IHE    WIHDCW    COHSTRAINTS 

The  time-ccnstrained  Traveling  Salesman  Problem  is  a 
variation  of  the  TSP  that  includes  time  window  constraints 
on  the  time  to  visit  some  of  the  cities.  The  hard  time- 
constrained  TSP  is  to  find  the  minimum  cost  tour  subject  to 
visiting  each   city   within   its   time    window. 

"For  the  time-constrained  TSP  model ,  we  define  a  contin- 
uous nonnegative  variable,  t^  ,  to  be  the  time  that  the 
salesman  visits  city  i.  Since  the  salesman  must  return  to 
city  1  (depot)  at  the  end  of  the  tour,  the  formulation 
includes  an  additional  variable,  tn+,  ,  the  total  time 
reguired  to   complete    the    tour- 

We  assume  that  a  complete,  symmetric,  nonnegative 
distance  matrix,  |Cij  J,  is  known  and  that  time  is  a  scalar 
transformation  cf  distance.  Thus  time  and  distance  may  be 
used   interchangeably. 

The  following  mathematical  formulation  of  the  TSP  with 
time  constraints   is   from    Baker  [2ef.    4]. 

MIN  t  -         t  (1.6) 

n+1  1 

S.T  t        -      t      >    c  i    =   2,3, ,n       (1.9) 

i  1  1i 

i    t      -    t      \    >    c  j    =    3,4,.. .,n      (1.10) 

2    <    i   <    j 

t      -      t        >    c  i    =    2,3, ,n      (1-11) 

n+1         i  1i 

t         >    0  i    =    1 ,2, ,n+1     (  1. 12) 

i 

1      <    t      <    u  i    =    2,3, ,n       (1-13) 

i  i  i 

where        t         =    the   time   that    the    salesman    visits   city  i 
i 

|x|       =    the   absolute   value   of    x 

c        =    the    shortest   time   reguired   to   travel    from 

city   i    to    city    j 
13 


1        =   the   lower  bound   on    the    time   window    for    the 
i 

salesman   to   visit  city   i 

by  assumption    all   1   >    0 

i 
u        =    the   upper   bound  on    the   time   window    for   the 
i 

salesman    to    visit  city   i 

u      >    1      ,    for    all   i 
i  i 

The  constraints  (1-9)  through  (1-12)  ensure  a  nonnega- 
tive  arrival  time  at  city  i,  t-  ,  be  obtained  for  each  city 
(node  2  through  node  u)  on  the  tour-  The  constraint  (1-9) 
guarantees  that  t;  ,  the  time  that  the  salesman  leaves  the 
node  1  (depot)  will  be  the  smallest  t  value.  The  absolute 
value  constraint  (1.10)  ensures  that  the  arrival  times  of 
any  two  city  differ  by  amount  of  time  sufficient  to  allow 
the  salesman  to  travel  between  the  two  cities.  The 
constraint  (1.11)  guarantees  that  t^  ,  the  time  the  salesman 
returns  to  the  depot,  will  be  the  largest  t  value-  The 
inegualities  (1.12)  and  (1.13)  are  nonnega  ti  vi  ty  and  the 
time  window  constraints  respectively. 

Unf ortuna tely.  Baker's  proposed  model  for  the  time- 
constrained  TSP  is  very  difficult  to  solve,  because 
constraint  (1.10)  is  nonconvex.  Therefore,  we  will  not  use 
this  formulation  directly  in  cur  program- 
Figure  1.1  illustrates  the  nonconvexity  of  constraints 
(1-10)  for  one  i,j  pair,  the  example  |  t  -  t  |  >  5.  The 
feasible  region  for   this  constraint   is  the   union  of   two 

disjoint   sets.    Taken   all   together,   constraints   (1.10) 

m 
define   2   disjoint   sets   sets    where   m   =   (n-1)  (n-2)/2, 

which  are  very  difficult  to  work  with. 

We   can   see   that   the   time-constrained   TSP   is   very 

different  from  ISP,  even  in  formulation. 


Mi 


10  ,- 
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Figure   1-1 


Example    of  Nonconvexity   of    (1.10) 
in   Two    Dimensions- 
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II.  HEDBISTIC  TSP  SOLOTION 

A.   OVERVIEW 

Many  heuristic  procedures  have  been  developed  for 
solving  TSP.  Our  purposes  in  this  Chapter  are  to  examine 
some  of  the  well-known  heuristics,  to  review  Stewart's  [  Hef - 
6]  recent  heuristic,  and  to  compare  these  approximate  tech- 
niques en  the  basis  of  efficiency  and  accuracy  on  a  small 
number    of   exanples. 

In    general,  heuristic    procedures      are   categorized      by 

three  broad  classes:  tour  construction  procedures,  tour 
improvement  procedures,  and  composite  procedures  [Ref-  7  ]. 
Tour  construction  procedures  start  with  a  single  node  and 
successively  add  nodes  till  a  tour  is  built.  Tour  improve- 
ment procedures  attempt  to  find  a  tetter  tour  given  an 
initial  tour.  Composite  procedures  construct  a  starting  tour 
from  one  of  the  tour  construction  procedures  and  then  try  to 
find  a  better  tour  using  one  cr  more  of  the  tour  improvement 
procedures. 

1-       Tour    Construction   Procedures 

There  are  many  methods  available  for  constructing  an 
initial  tour.  Procedures  which  have  been  generally  used  are 
given    below. 

a.       Nearest   Neighbor     (   Rosenkrantz    et   al.    [Ref.    3]    ) 

Step    1.    Start    with    any   node   as   the   beginning   of 

a      subtour. 
Step   2.    Find    the   node   closest    to   the    last    node 

added   to    the    subtour.    Add    this    node    to 

the   current    subtour. 
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Step  3.    Repeat   step    2   until    all   nodes    are   contained 
in    the   tour.    Then,     join   the   first   and    last 
node. 

b.  Clarke   and      Wright  Savings       (   Clarke      and   Wright 
[fief.    9]   ) 

Step    1.    Select   any   node   as    the   central    depot 

which   we   denote   as    node    1. 

Step   2.    Compute    savings  s        =   c        ♦      c  c 

ij     1i      Ij      ij 

for  i,  j  =  2,3,...,n.  i*j 
Step  3.  Order  the  savings  from  largest  to  smallest. 
Step  4.  Starting  with  the  largest  savings  on  the 
list,  subtours  are  assembled  such  that  the 
next  node  added  has  the  largest  remaining 
savings  -  provided  that  a  constraint  is  not 
violated.  Once  a  pair  of  nodes  i  and  j  have 
been  linked,  they  remained  linked  . 
Repeat  until  all  nodes  have  been  assigned. 

Here,  tne  guantity  s.  is  the  amount  of  travel 
saved  if  node  j  is  visited  directly  after  i,  as  opposed  to 
having  separate  trips  from  the  depot  to  nodes  i  and  j. 
Figure  2.1  demonstrates  the  procedure  for  two  nodes  i  and  j. 

c.  Insertion  procedures  (  Rosenkrantz  et  al-   [Ref. 
8]  ) 

Ac  insertion  algorithm  constructs  a  feasible 
tcur  by  successively  adding  one  node  to  an  existing  subtour. 
This  procedure  takes  a  subtour  of  k  nodes  at  iteration  k  and 
attempts  to  determine  which  node  not  in  the  subtour  should 
join  the  subtour  next  (the  Selection  step).  Then  it  deter- 
mines where  in  the  subtour   it  should  be  inserted  {Insertion 
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Figure   2.1        Concept    of   the   Clarke  -   Bright 
Savings    Heuristic. 


step).       Stewart      [fief.      6]      presented    the      following   general 
algorithmic    structure- 
Step    1.     (Initial      Suttour) 

Obtain  a   IS P   tour   for   a    subset    of    the    nodes 
N 'C  N    in    G- 
Step   2.     (Selection    Step) 

Find    a   node    k €  N-N '    tc   be    added    to    the 
existing    suttour. 
Step   3.     (Insertion    Step) 

Choose    an   arc(i,j)     in    the   subtcur    on    N*. 
Insert   node    k   between    i   and    j    and      add 
k    to   N' . 
Step   4.    If    N   =    N«,    then    stop- 

(We    have   a    HamiltoniaE  cycle). 
Otherwise,    return    to    step    2. 
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There  are  many  variations  on  this  algorithmic 
structure  depending  on  the  procedures  chosen  for  executing 
steps    1  r2   and  3  . 

Wiorkowski  and  McElvain  [Ref.  10],  Or  [fief.  11], 
Stewart  [Ref-  12]  and  Norback  and  Love  [fief.  13]  all  present 
insertion  algorithms  that  use  the  convex  hull  of  the  set  of 
nodes  N  for  the  initial  subset  N1.  Nemhauser  and  Hardgrave 
[fief.  14]  have  shown  that  there  exists  an  optimal  tour  for 
every  Euclidean  TSP  in  which  the  relative  order  of  the  nodes 
on  the  boundary  of  the  convex  hull  is  preserved.  This  means 
that  the  optimal  tour  visits  codes  on  the  boundary  of  the 
convex  hull  in  the  same  order  as  if  the  boundary  itself  were 
followed. 

Further  justification  for  the  use  of  the  convex 
hull  for  the  initial  subtour  is  shown  empirically  by 
Stewart's  [Ref.  6]  computational  experiment.  He  compared 
several  insertion  heuristics  both  with  and  without  the 
convex  hull  as  the  starting  solution.  The  results  show  that 
all  the  insertion  algorithms  are  improved  by  the  use  of  the 
convex  hull.  Some  are  improved  substantially,  others  only 
moderately. 

Many  criteria  have  been  suggested  for  the  selec- 
tion  of   the    node   to   be   inserted   in    an    insertion    procedure. 

( 1)  Nearest    Neighbor    (    Rqsenkrantz    et    al-    f  Ref. 
8]   ).        Choose    the    node  k    that   is   nearest   a    node    in   the 
current    tour.    I.    e-    ,    find      k    =   argmin      c  s.t-    j€N-N«, 

ieH..  j       ij 

(2)  Cheapest      Insertion       (   fiosenkrantz      et      al. 

[Jtef-       8]      ) •         Choose    the      node    k      that    may    be      inserted   at 

minimal     increased   cost.      I.e.,      find 

k  =   argnin      c        +c        -c  s.t.    me  N-N1  ,  i,  j  €  N* . 

m  im  mj  ij 
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(3)   Farthest   Insertion   (  Bosenkrantz   et   al. 
r  Ref.  8])-    Choose  the  node  k  that  is  farthest  from  a  node 
current  subtour.   I.e.,  find   k  =  argmax  c     s.t.  j€N-N*, 

ie«.  3      ij 

(*4)  Arbitrary  Insertion  (  Rosenkrantz  et  al. 
r  Ref.  8]  ).    Choose  node  k  randomly  from  among  N-N*. 

(5)  Ratio  Insertion   (  Stewart  [Ref.    J2  ]  ). 

Choose  the  node  k  such  that  the  proportional  increase  in 

ccst  is  minimal.  I.e.,   find  k  =   argmin  (c   ♦  c    )  /  c 

m  im  mi  ij 

s.t-    m   6   N-N» ,        i, j  €  N«  - 

(6)  Perpendicular  Distance  (  Wiorkowski  and 
McElvain  [Ref-  JO]  ) . -  Choose  the  node  k  that  is  closest  to 
an   arc    in   the  current    subtour. 

(7)  Ratio   Times      Distance    (    Or   [Ref.         11 1    ). 
Choose        the    node      x      such   that      the      product      of    ratio      and 
distance  is      minimized.      I.    e.    ,    find 

k   =    argmin    ((c        +   c      )/c      )x(c        +c        -c      ) 

m  im  mj  ij  im  mj  ij 

s.t.  IE    €    N-N'     ,         i,  j      €     N  ». 

(8)  Greatest   Angle         (    Norback   and       Love      [fief- 
131    )-         Choose    the   node    k  and     arc    if j      such    that    the    angle 
formed    by    the   two   arcs    (i,k)     and    (k,j)    is   a   maximum.      I.e., 
find   k    =   argmax      angle{   arc(i,m),    arc{m,j)    }       s.t.    m€N-N', 

ID 

i,j    €    N«    . 

The    insertion   criteria   that    have   been    used 
fall   into   two    categories.    [Ref.    6] 

1.       Cheapest   Insertion 

Insert    the    node   kfc  N-N*    between    those    two 

connected      ncdes   i,j    €    N'      that      minimize   the      guantity 

c        +    c        -    c 
ik  kj  ij 
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2.    Identical   Insertion    and   Selection 

Do      selection   and      insertion      in   the      same 
step. 

2-       Tour    Improvement    Procedures 

The  best  known  procedures  of  this  type  for  the  TSP 
are  the  branch  exchange  heuristics  [Ref.  7].  These  branch 
exchange  heuristics   work    as    follows. 

Step   1.    Find   an    initial    tour.    This   tour   may    be   chosen 
randomly    from   the   set   of    all      possible    tours, 
cr   it   may      be   generated      by    one   of      the   tour 
building    procedures    above. 

Step  2.    Improve      the    tour      using   one      of    the      branch 
exchange    heuristics. 

Step   3.    Continue    step   2,      until   no      additional 
improvement    can    be    made. 

For  a  given  kr  we  define  a  k-change  of  a  tour  as 
consisting  of  the  deletion  of  k  branches  in  a  tour  and  their 
replacement  by  k  other  branches  to  form  a  new  tour.  A  toar 
is  k-opt  if  it  is  not  possible  to  improve  the  tour  via  a 
k-change.  In  general  the  larger  the  value  of  k,  the  more 
likely  it  is  that  a  k-opt  solution  will  be  optimal. 
Unfortunately,  the  number  of  operations  necessary  to  test 
all  k  exchange  is  proportional  to  n  r  where  n  is  the  number 
of  nodes  in  the  TSP.  Due  to  this  complexity,  values  of  k  =  2 
and  k  =  3  are  most  commonly  used  [  Ref .  7].  The  2-opt  and 
3-opt   heuristics   were      introduced   by   Lin    [Ref-  15]    and   tiie 

k-opt  procedure,  for  k> 3  was  presented  by  Lin  and  Kernighan 
[ Kef.    1 6]    . 

Or  [Ref.  11  ]  has  designed  a  modified  3-opt  that 
considers  only  a  small  percentage  of  3-branch  exchanges. 
This      modified      3-opt      called   Oropt      by      Stewart      [Ref-  6] 
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considers  onlj  those  branch  exchanges  which  are  composed  of 
a  string  of  one,  two,  or  three  adjacent  nodes  being  inserted 
between  two  ether  nodes  in  the  current  tour.  3y  limiting  the 
number  of  exchanges  that  are  considered  in  this  way,  Oropt 
reguires   many   fewer   calculations    than    a   full   3-opt. 

Stewart  [Eef.  6]  made  an  experiment  of  the  convex 
hull,  cheapest  angle  insertion  algorithm  (CCA)  which  will  be 
discussed  in  the  next  section  as  a  stand-alone  algorithm  and 
with  each  of  the  three  post- processors.  The  algorithms  are 
designated  CCA,  CCA2,  CCAO,  and  CCA3  for  the  convex  hull 
cheapest  insertion  stand-alone,  with  2-opt,  with  Oropt  and 
with  3-cpt  respectively.  He  drew  two  conclusions  from  his 
computational  results.  First,  the  3-opt  reguires  substan- 
tially more  time  than  either  the  2-opt  or  tne  Oropt. 
Second,  the  2-opt  is  dominated  by  the  Oropt  and  the  3-opt  in 
guality    of    solution. 

In        computation        time,         Oropt        only        looks        at 

2 
approximately      3n        of    the      n    possible      3-opt      exchanges      on 

eaca   pass.       There   are      n    ways      to      select    tne    first      branch, 

times    3    ways    to    select   the      second    branch,         and    n-2   ways   to 

select    the    third    branch. 

This  accounts  for  the  fairly  close  times  for  the 
2-opt   and   Oroft.  The      guality    of    CCAO      solutions    dominate 

CCA2  solutions.  On  the  other  hand,  there  is  little  or  no 
difference  between  the  Oropt  ani  3-opt  in  terms  of  solution 
guality - 

Stewart's  main  conclusion  from  the  above  experiment 
is  that  the  Orcpt  performs  as  well  as  a  3-opt  in  a  small 
percentage  of  the  computer  time  reguired  by  a  3-opt,  and  it 
should  be  preferred  tc  both  the  2-opt  and  the  3-opt  for 
Euclidean   TSP's. 
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3.      Composite    Procedure 

The  basic  composite  procedure  is  a  combination  of 
the  tour  construction  and  branch  exchange  procedures.  It  is 
obtained  by  appending  a  branch  exchange  procedure  to  the 
tour  construction  algorithm  as  a  post- processor.  The  proce- 
dure  can  be    stated   as   follows  [  Ref .     17]. 

Step    1-    Obtain  an  initial   tour    using    one    of    the 

tour  construction  procedures. 
Step  2.    Apply   a  branch   exchange   procedure   to   the 

solution    produced  by   the    step    1. 

Stop   wnen    no   further   improvement    can    be   made. 

The  composite  procedure  is  relatively  fast  computa- 
tionally and    gives    good   results  [Ref.     18]. 

E-      CCAC 

1  •       Algorithm 

The  GCAO  algorithm  designed  by  Stewart  [fief.  6]  uses 
the  convex  hull  of  the  nodes  in  N  for  its  initial  subtour. 
Then  it  inserts  the  nodes  not  currently  in  the  subtour  where 
they  may  be  inserted  most  cheaply  (the  Cheapest  Insertion 
criterion).  It  selects  the  node  k  tc  be  inserted  at  each 
iteration  according  to  how  large  an  angle  is  formed  by  the 
two  arcs  that  must  be  added  to  the  current  subtour 
(Selection  criterion)  in  order  to  insert  k.  Finally  it  uses 
an  Oropt  to  make  local  improvement  on  the  tour  constructed 
in  the  first  stage.  CCAO  means  Convex  Hull,  Cheapest 
Insertion,    Angle   Selection,    Oropt. 

Algorithm    :    CCAO 

Input    :    Number    of    nodes,      x   and    y   co-ordinates   of    all 

nodes. 
Output:    Ordered   list    of    tour,    total   cost. 
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Step    1    .     (Initial    Subtour) 

Find   the    convex   hull    of   the   set    of   nodes   N. 

Call    the      set   of    nodes   on    the   boundary      N1. 

Let   the    initial   subtour   be  the      nodes    of    N* 

in   the      same      order    as     they      appear   on   the 

convex   hull. 
Step    2    .     (Cheapest  Insertion) 

For   each   node   m€N-N',    find 

(i    ,  j   )     =   argmin      c  ♦   c        -    c 

m      m  i,j  im  mj  ij 

s.t.   i,  j  €    N* ,   ir  j  :  connected. 

Step  3  .  (Greatest  Angle  Selection) 

For  the  next   insertion,  select   the   node 

that   maximizes  the  angle  between  the  arcs 

(i  ,  m)  and  (m,j  )   over  all  m  6  N-N'. 
m  m 

I.e,    find   k    =   argmax    angle {    (i    ,m),(a,j    )    } 

m  mm 

s.t.         m€  N-N». 

Insert      k         tetweeh    i        and    j         and   add 

k       k 
k   to  N' . 

step  4  :   If  N«  =  N,  go  to  step  5- 

Otherwise  return  to  step  2. 

step  5  ;   Apply  an  Oropt  to   the  current  tour.  Stop 

when  no  further  improvements  can  be  found. 

End  of  algorithm  CCAO 

2 .   Exa  mple 

Figures  2.2  -  2-  4  illustrate  the  above  algorithm  on 
the  TSP  defined  as  test  problem  [ 1 ]  in  Appendix  A.  First 
the  convex  hull  is  generated  for  an  initial  starting 
subtour.  This  subtour  consists  of  nodes  2,13,12,14,5,15,7,4.. 
A  solid  line  marks  the  boundary  of  the  convex  hull  in  Figure 
2.2. 
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Figure  2.2    Initial  Sabtoar  and  Insertion. 

In  step  2,   each  of   the  interior  nodes  (1,3,6,6,9,10,11,16) 
is  associated  witn  a  pair  of   connected  nodes  on  the  initial 
suntour  (the  dashed   lines  in  Figure  2.2).   In   step  3,   the 
dashed  lines  that  form  the   greatest  angle  (closest  to  180°) 
identify  the  node  to  be  inserted  (node  10   in  this  example). 

Figure  2.3  shows  the  problem  after  the  first  three 
insertions (  ncde  10,  node  1  and  node  8  in  that  order). 
Notice  that   some  nodes   not  in   the  subtour   are  associated 
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Figure    2.3         Internediate   Subtour   and   Insertions. 

with  new  node  pairs.  Figure  2.4  shows  the  final  tour  for 
stage   ore.  This    tour      is      now    passed      to      an    Oropt      post- 

processor. In  this  case  the  tour  from  stage  one  appears 
from  inspection  to  be  optiial,  and  Oropt  will  find  no 
improve  irent . 
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Figure   2-4        Final   Tour   of  CCAO. 

3.      Computational    Results 

In  addition  to  CCAO,  CCCO  (Convex,  Cheapest, 
Cheapest,  Oropt)  has  been  coded  for  the  purpose  of  compar- 
ison. The  only  difference  fcetveen  CCAO  and  CCCO  is  that 
CCCO  uses  the  cheapest  selection  criterion  instead  of 
greatest   angle    of    CCAO. 

We  used  SedgevicJc's  [  Eef .  19]  package  wrapping  algo- 
rithm     for      finding      the        convex      hull      (initial      subtour) - 
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Starting  with  some  point  (called  the  anchor)  that  is  guaran- 
teed to  be  on  the  convex  hull  (say  the  one  with  the  smallest 
y  co-ordinate),  take  a  horizontal  ray  in  the  positive  direc- 
tion and  sweep  it* upward  until  hitting  another  point-  This 
pcint  i s  on  the  hull.  Then  start  at  that  point  and  continue 
sweeping  until  hitting  another  point,  etc.  The  package  is 
completely  wrapped  when  the  first  point  is  included  again. 
The  following  algorithm  finds  the  convex  hull  of  an  array 
1(1,..., a)  of  nodes,  the  node  L(n+1)  is  used  as  a  sentinel, 
that  is,  a  copy  of  the  first  node  which  is  used  to  signal 
completion  of  the  procedure.  The  variable  NH  is  maintained 
as  the  number  of  nodes  so  far  included  on  the  hull- 
Algorithm  z    Package  Wrapping 

Input  :  Number  of  nodes,  x  and  y  co-ordinates  01  ail 
nodes. 

Output:  Ordered  list  of  convex   hull  and   number  of 
nodes  included  on  the  convex  hull- 


Step  1  .  (Initialization) 

find  and   duplicate  anchor.    I.e.,   find 

NHIN  =   argmin   y    s.  t.  i  €    K   and    set 

i 

NH    =    0,       l(n  +  1)    =    L  (NMIN)  . 

Step    2    ;     (Swap    nodes    NH   and    NMIN). 

Put      last      node      found    into      the      hull      by 

exchanging   it    with    the    NHth    node. 

NH    =    NH    ♦    1. 

TEMP    =    L  (NH)  . 

L(NH)     =    L  (NMIN)  . 

I(NMIN)    =    TEMP. 
Step    3    :     (Compute    angle) 

Compute   the    angle    from      tne    horizontal    made 
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by    the   line    between    L(NH)    and   each   of      the 
nodes   not   yet   included    on   the   hull- 
Step    4    :     (Find    next  hull   node) 

Find   the   node    whose    angle  is   smallest   amcng 
those    with      angles      bigger    than    tne   current 
value  of    the      'sweep'    angle    (the   angle    from 
horizontal      to   the      line        between      L(NH-1) 
and   L(NH)  )  - 
Step    4    :    Stop      when   the  first      point   is      encountered 
again.         I.e.,      L(n+1)     =   L(NHIN)     . 
Otherwise,   go    to   step    2. 
End  of  algorithm  Package  Wrapping 

We  used  Sedgewick's  Pseudo  Angle  for  finding  the 
smallest  angle  in  step  3,  which  is  coded  as  the  'THETA' 
function.  This  function  returns  a  real  number  between  0.0  to 
4.0  that  is  net  the  angle  made  ny  L1  and  L2  with  the  hori- 
zontal but  which  has  the  same  order  properties  as  the  true 
angle.  If  dx  and  dy  are  the  delta  x  and  y  distances  from 
some  node  to  the  anchor  node,  then  the  angle  needed  in  this 
algorithm  is  arctangent  dy/dx.  However,  the  arctangent 
function  is  likely  to  be  slow  and  it  leads  to  at  least  two 
annoying  extra  conditions  to  compute  :  whether  dx  is  zero, 
and    which   guadrant    the    point    is   in. 

In  this  algorithm  we  only  need  to  be  able  tc  compare 
angles,  not  measure  them.  Thus  it  makes  sense  to  use  a  func- 
tion that  is  much  easier  to  compute  than  the  true  angle  but 
has  the  same  crdering  properties  as  the  true  angle-  A  good 
candidate  for  such  a  function  is  siaply  dy  /  (dy  +  dx)  . 
Testing  for  exceptional  conditions  are  still  necessary,  but 
simple  r. 

Function   THETA    (   Pseudo    Angle  ) 

Input      :    dx,dy    (delta    x  and    y   distances   from    some 
node    to   the    anchor    node) - 


29 


Cutput   :    Pseudo      angle    made   by    11    and  L2    with   the 
horizontal  line. 

begin 

dx  =  x  (L2)  -  x(L1)  :  ax  =  abs(ax)  : 
dy  =  y(12)  -  y(L1)  :  ay  =  abs  (ay)  : 
if    (    dx=0   )    and    (    dy=0    )       then    t  =   0.0 

else    t   =   dy  /    (ax   ♦   ay    ) 
if   dx    <    0      then      t   =    2.0    -   t 

else   if   dy    <  0    then    t   =   4-0    +    t 
end 

End   of    function   THETA 

Figure  2.2  shows  how  the  hull  is  discovered  in  this 
way.  We  used  Sedgewick's  Pseudo  Angle  for  finding  the 
greatest  angle    selection    point. 

The  data  for  our  test  problems  is  given  in  the 
Appendix.  The  computational  results  are  summarized  in  Table 
I.  As  can  be  seen  in  Table  I,  CCAO  is  taster  than  CCCO  on 
the  small-scaled  test  problems  (below  30  nodes  ) ,  nut  CCCO 
is  faster  than  CCAO  on  the  moderately  large  sized  problems 
(over  5  0  nodes) .  Generally,  the  accuracy  is  almost  iden- 
tical  in  both   cases. 

Stewart  [Ref.  6]  showed  that  in  a  large  scaled 
problem,  the  CCAO  algorithm  outperforms  any  other  insertion 
and  selection  algorithms.  Thus,  we  are  highly  motivated  to 
use  a  modification  of  the  CCAO  algorithm  for  solving  the 
time-window    constrained    TSP- 
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TABLE  I 
CCHPOTATIOBAL  BESDLTS  OF  CCCO,  CCAO 


I 

CCCO 

I 

CCAO                  J 

Problem 

Number 

Number 

of 
Nodes   n 

Best 

Known 

Solution 

% 
Over 
Best 

CPU 
Time 
(sec) 

or 

Over 
Best 

CPU 

Time 

(sec) 

[1] 

16 

66.603  9 

0.00 

0.0133 

0.00 

0.0066 

[2] 

22 

469.0288 

0.00 

0.0233 

0.00 

0.0100 

[3] 

22 

278.4371 

0.00 

0.0166 

0.00 

0.0066 

[5] 

51 

429.7000 

2.72 

0.  1897 

3.  94 

0.2562 

[6] 

76 

552.9000 

1-  64 

0.5857 

1.54 

0.6889 

.    i 

*  CPU  times  in  seconds  on  IBM  3033- 
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III.     1JE    ISP    WITH    HAR D    TIHE    WINDOW    CONSTRAINTS 

A.       INTBODUCTION 

The  first  time-constrained  TSP  we  consider  is  the  case 
in  which  late  arrivals  are  not  allowed,  and  early  arrivals 
must  wait  for  the  opening  of  the  time  window  before  they  can 
begin  to  service  a  customer-  This  is  called  the  hard  time 
window    case    and    it    is    illustrated   in   Figure    3.1. 


waiting 
—  — >j  time  l<  — 

salesman  lower 

arrives  bound 


cannot   be 
violated 


J__ 

upper 
bound 


time   window 
for  city  i 


.     Figure    3- 1         Diagran   for   Hard    Tiie  Window    Case. 

The  hard  time  window  case  corresponds  to  military  opera- 
tions and  to  some  civilian  distribution  problems.  Meeting  a 
deadline  is  considered  a  critical  factor  in  this  case.  The 
soft    time   window   case    will  be  .discussed   in    the    nextChapter. 

Consider  a  graph  G  =  {N,A}  composed  of  a  set  of  nodes  N 
and  a  set  of  arcs  A  connecting  these  nodes.  We  now  define 
some  notation  to  be  used  throughout  cur  discussion  of  the 
time-window-ccnstrained   TSP. 
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1  =  Lower  bound  on  the  time  window  at  node  i 
(early  allowable  arrival  time  at  city  i)  . 

u    =  Upper  bound  on  the  time  window  at  node  i 

(latest  allowable  arrival  time  at  city  i)  . 

d    =  Time  required  to  spend  at  node  i. 
i 

(service  time  at  city  i) - 

SPEED  =   Constant  speed   at   which  the   vehicle 

travels. 

dist  =   Distance  from   i   to   j. 

ij 

c  =   Iravel    time  from   i    to    j. 

ij 

Note    :    c        =      dist  /    SPEED. 

ij  ij 

We   use   c  and     c(i, j)    interchangeably, 

ij 

depot   =    Depot (home)     node. 

I  =     (    L(1)  ,L(2)  ,...,L(n)    ). 

=  A  tour  with  n  stops  visited  in  the  order 

1(1), 1(2)  ,.-.,L  (n)  . 

ARRV1      =      Arrival   time   at   city    i. 

i 
MAI T      =    Waiting      time    at    node    i      for    the    hard   time 
i 

window. 

He    also    use   l(i),    u(i),    d(i)r.  ARRVT(i),    WAIT(i)     and 

1    r    u      r    d       t    ARRVT       ,    WAIT      interchangeably. 
i        i  i  i  i 


B.       HEURISTIC    SOLUTION    TECHHIQUES    FOB    BARD    TIME    WINDOWS 

1.       Nearjst    Neighbor 

The  following  is  a  Nearest  Neighbor  heuristic 
similar  to  the  one  used  in  the  unconstrained  TSP.  At  each 
iteration  we  add  a  new  node  tc  the  end  of  the  subtour.  It 
is   the    first    node    that  can   be      visited    from    the    last   node    of 
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the    sub  tour ,    taking  into    account    any   waiting   time   that    might 
te   necessary    due   to  the   lower   time    window  bounds. 

Algorithm    :    Nearest   Neighbor 

Input       :    Number    of    nodes,    x  and    y   co-ordinates  of    all 

nodes,    time   windows    for    all   nodes- 
Output    ;    Ordered    list    of   tour,    total    travel   time. 

Step    1    .     (Initialization) 

Start    at   the    depot. 

Let    i=depot,       N«    =    {i} . 
Step    2    .    Compute   ARRVT    for    all  nodes      k€  N-N»    if    k   can 

be   visited    directly   after   i    : 

AERVT    =    max    {    AERVT    ,    1   }    +    d    +    C 

It  i   i     i   ik 

Step  3  .  If  ARRVT   >  u  ,   then  stop  ('  no   feasible 

k     k 

solution1 )  . 

Step  4  .  If  AERVT  <  1  ,  then   cost  =   1   . 

k        k  k  k 

Otherwise,      cost  =    ARRVT 

k  k 

Step    5    -     (   Nearest    Neighbor   Selection) 

Choose    the   node   k€  N-N'    such   that      cost      is 

k 
a    minimum.      I.e,      find 

k   =   argmin     cost        s.t.    j  6N-N'. 

J 

Step  6  .  (Insertion) 

Insert  k   after  i,  add  k   to  subtour  N',  and 

let  i  =  k- 
Step  7  „  If  N'=Nr  go  to  next  step. 

Otnerwise,  return  to  step  2. 

Step  8  »  Compute  total  travel  time,  then  stop. 

Total  travel  time  =  max  {  ARRVT  ,  1   } 

k         k 

♦    d    +   c 

k        k, depot 

End   of   algorithm   Nearest   Neighbor 
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This  solution  was  constructed  by  starting  at  the 
depot  and  moving  to  the  nearest  neighboring  customer  that 
has  not  yet  teen  visited  as  long  as  the  upper  bound  level 
was  not  violated.  This  heuristic  may  fail  to  solve  the 
problem . 

2.   SCCO 

This  algorithm  is  designed  for  the  case  when  some  of 
the  nodes  do  not  have  time  windows.  We  call  these  nodes  " 
time   free  ". 

SCCO  is  similar  to  the  cheapest  selection,  cheapest 
insertion  method  for  the  unconstrained  TSP,  except  that  the 
nodes  with  time  windows  are  treated  differently  from  the 
time  free  nodes.  The  nodes  with  windows  are  inserted  in 
order   of  increasing    upper    time   window    tound. 

The  time  free  nodes  are  inserted  between  these  nodes 
by  cheapest  selection  and  cheapest  insertion,  for  as  long  as 
the    upper      bound   time      window    will   allow.  In    the      end,      a 

Modified  oropt   is   used   to   improve   the   solution. 

There  is  one  possible  difficulty  with  this  approach. 
It  may  become  impossible  to  reach  some  of  the  tiae- 
ccnstrained  nodes  before  their  upper  bound.  In  this  case, 
we  must  delete  some  node  (s)  from  the  subtour.  Whenever  we 
see  an  upper  oound  that  cannot  be  satisfied,  we  select  a 
node    to    delete    ty    the   following  criteria. 

The  first  criterion  is  the  width  of  the  time  window. 
Hence,      time-free      nodes   are      considered   first.  Then,       if 

several  nodes  in  the  subtour  are  tied  for  the  widest  time 
window,  we  select  for  deletion  the  node  that  results  in  the 
greatest     time      saved.  The      algorithm      is      summarized      as 

fellows . 

Algorithm    :    Successive  Cheapest  Cheapest  Oropt    (SCCO) 

Input       :    cumber    of    nodes,    x  and    y   co-ordinates  of    all 
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nodes,  time  windows  for  all  nodes. 
Output  :  ordered  list  of  tour,  total  travel  time- 
Step  1  .  (Initialization) 

Start  at  the  depot. 

let   i=depot,  N«  =  {i}. 

Step    2    .    Set    k    =   argmin    u        s.t.         J6N-N1. 

J 

If   k   is   time    free    node,    then    set    k  =    depot. 

Step    3    .    Calculate    ARRVT  . 

k 

ARRVT  =  max  {  ARRVT  ,  1  }  +  d  +  c 

k  i    i     i    ik 

Step  4  .  If  ARRVT   <  u   ,   then  go  to  step  5. 

k     k 

Otherwise,  select  time  free  node  m  e  N1  which 
results  in  the  greatest  time  saved  for  dele- 
tion. Delete  node  m  from  N*,  go  to  step  3. 

Step  5  .  Add   node  k  to  the  subtour  N*. 

Step  6  .  Insert  time  free  node  j  €  N-N1  between  nodes  i 

and  k   by   cheapest  insertion   and    cheapest 

selection  (same  as  CCCO)  until  ARRVT  does  not 

k 

exceed  u 

k 

If  ARRVT   <  1  ,  then   set  ARRVT  =  1  . 
k     k  k    k 

Step  7  -  If  N'=N,  then  go  to  next  step. 

Otherwise,  let  i  =  k  and  go  to  step  2. 

Step  8  .  Apply  the  Modified  Oropt   procedure  to  the 
current  tour.  Stop  when  no  further  improve- 
ments can  be  found. 

End  of  algorithm  SCCO 

"Successive"  means   select  the  node   successively  by 
the  smallest  upper   bound.    In  the  SCCO   algorithm,   if  the 
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salesman  arrives  before  the  lower  bound  of  the  time  window, 
adding  waiting  time,  we  set  the  arrival  time  equal  to  the 
lower   bound. 

The  Modified  Oropt  procedure  for- improving  the  solu- 
tion is  described  below.  This  procedure  consider  only  those 
exchanges  that  would  result  in  a  node  being  inserted  between 
two   other  nodes   in   the  current  tour. 


X 


9        10       11       12    j 


Figure    3-2        Curreat  Tour  before   Hodified-Oropt- 

Figures  3-2  and  3.3  are  helpful  to  understand  how 
the  procedure  works.  In  both  figures,  i,j,k,l,  and  m  are  the 
nodes    in  the    current    tour.       Nodes      1   and   m    are    considered   to 
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be  adjacent  tc  node  Jc.  A  test  is  then  conducted  to  deter- 
mine if  node  k  can  be  located  between  tvo  other  nodes,  such 
as  i  and  j,  sc  that  it  results  in  reduced  total  travel  time- 
If  it  can,  we  make  the  appropriate  arc  exchanges,  then 
update    the    total  cost   and   route  orders. 
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Figure   3*3        Iaproved    Tcur   after   Modif ied-Oropt. 


In    this   example,      the    three    arcs    (i,j) ,  ^k,l) ,      and 

<k,m)    are   removed   and    replaced   by    (i#k),       (j#k),  and    (l,m). 

When   no    further      exchanges   improve    the    solution,  the    algo- 
rithm  terminates. 


38 


Figure  3.4        Subtour  in  SCCO   Procedure. 

Figures  3.  4  -  3.6  illustrate  the  SCCO  algorithm  for 
the  TSP  with  hard  time  windows  given  in  Appendix  F  as  in 
test  problem  [1].  In  this  problem  10  of  16  nodes  have  time 
windows.    The    ether   6    nodes  are   time    free. 

First,  the  subtour  starts  at  the  depot  (node  16)  and 
we  insert  the  node  with  the  smallest  upper  bound  (node  12). 
We  examine  all  nodes  which  could  be  inserted  between  16  and 
12   as   long      as   the    upper    bound      on    node    12    is      observed.      In 
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Figure  3.5   Intermediate  Subtoar  in  SCCO  Procedure. 

this  case  there  is  no  such  node.  Then  we  select  the  next 
smallest  upper  bound  (node  1 4) ,  add  it  to  the  tour,  look  for 
nodes  to  insert  before  it,  and  continue  in  this  manner.  Now 
we  have  formed  the  partial  tour  16,  12,  14,  11,  6,  3  as 
shown  in  Figure  3.4. 

As  shewn  in  Figure  3.5,  we  can  insert  3  time  free 
nodes  between  node  6  and  node  3.  These  insertions  are  made 
according  to   the  cheapest  insertion  and   cheapest  selection 
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Figure   3.6        Final  Subtour   for   SCCO. 

criteria.         He      do   not   make      any   further      insertions   because 
they    would    cause   a    time    window   violation   at   node    3. 

Figure  3.6  shows  the  final  tour  for  the  SCCO 
heuristic.  This  tour  is  passed  to  a  Modified  Oropt,  but  in 
this   case   it    will   find   no  improvement. 
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3.       SCAO 

This  heuristic  is  identical  tc  SCCO  except  for  the 
use  of  the  greatest  angle  selection  criterion  for  the  time- 
free   nodes,    instead  of   cheapest  selection. 

Algorithm    :    Successive  Cheapest   Angle  Oropt    (SCAO) 

Input       :    Numier    of    nodes,    x   and    y   co-ordinates  of    all 

nodes,    time  windows   for    all   nodes. 
Output    :    Ordered    list    of   tour,    total   travel   time. 

Step    1    .     (Initialization) 

Start   at   the    depot. 

Let    i=depot,     N*    =     {i}. 

Step    2    .    Set    k    =   argmin      u         s.t.         j    €    N-N'. 

D 

If   k    is   time    free   node,    set    k   =    depot. 

Step   3    .    Calculate    ARRVT   . 

k 

ARRVT    =    max    {    ARRVT    ,    1    }    +    d    +    c 

k  .i        i  i        ik 

Step    4    .    If    AivRVT      <   u       ,      then    go   to  step    5. 

k     k 

Otherwise,  select  time  free  node  m£N'  which 
results  in  the  greatest  time  saved  for  dele- 
tion. Delete  node  m  from  N',  go  to  step  3. 

Step  5  .  Add  node  k  to  the  subtour  N'. 

Step  6  .  Insert  time  free  node  j  €  N-N'  between  nodes  i 
and  k  by  cheapest  insertion  and  greatest 
angle  selection  (same  as  CCAO)   until   ARRVT 

K 

does   not    exceed      u 

k 

If    ARRVT      <    1    ,    then      set    ARRVT    =    I    . 
k  k  k        k 

Step    7    .Leti   =    k. 

If    N'=N,    then    go   to    next    step. 
Otherwise,    go    to    step    2- 
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Step    8    -    Apply    the    Modified      Oropt    procedure    to   the 
current  tour.     Stop   when    no    further    improve- 
ment  can   be  found. 

End  of  algorithm  SCAO 

This  algorithm  is  same  as  SCCO  except  greatest  angle 
selection   is    used   instead   of   cheapest    selection,    as   in    SCCO. 

4.       SLACK 

This  heuristic  was  suggested  by  Professor  Rosenthal. 
It  is  designed  for  the  case  when  the  widths  between  the 
upper  and  lower  bounds  of  the  time  windows  are  relatively 
large. 

In  this  heuristic,  the  SLACK  is  the  most  important 
concept.  The  SLACK (i)  can  be  defined  as  the  maximum  amount 
of  time  by  wnich  arrival  at  node  i  can  be  delayed 
without  causing  an  upper  bound  to  ne  violated  for  a  node 
currently   on    the   tour. 

The  SLACK  function  can  be  defined  as  a  recursive 
function   as    fcllcws. 

SLACK(L(i))    =    min     {    u(L(i))       -    AEfiVT(L(i))     , 

SLACK  (L(i*1))     +    WAlT(L(i))     } 

where 

WAIT(L(i))    =      niax     {0,    1(1  (i))     -    AKEVT  {L  (i)  )     } 

The  first  element  of  this  recursive  function  is  the 
difference  between  the  upper  bound  and  arrival  time  at  node 
L(i).  The  second  one  is  the  sum  of  next  node's  SLACK  and 
waiting      time      ex    node      L  (i)  .  The      minimum    of      these      two 

elements  is       a    possible      delay    time    of      the   arrival      time   at 
ncde  L ( i)  without   violating   tne      upper    bound    of      all    nodes 

after   L  (i)    in  the   current    tour. 
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The  advantage  of  this  recursive  function  is  that  it 
is  easy  to  calculate  a  possible  delay  time  without  calcu- 
lating new  arrival  times  for  all  nodes  after  L  (i) -  Ihe  algo- 
rithm  is  summarized   as  follows. 

Algorithm    :    SLACK 

Input    :    Number    of   nodes,    x   and    y   co-ordinates   of    all 

nodes,    time    windows   for    all   nodes- 
Output    :    Ordered   list    of    tour,    total    travel    time. 

Step    1    .     (Initialization) 

Start   at    the    depot-    Let    N»    =    {depot}. 

Step    2    -    Sort    the    upper   time    windows. 

U   =    (    u   ,u    ,.  .  .,u    ) 
1      2  n 

s.t.         u   <  u  <  .  ..<u 
1      2  n 

Step    3    .    Set    k    =   argmin    u        s-t.    i€  N-N1. 

i 

Step  4  .  Find  a  node  L(ISTAR)   after  waich   node  k 
can  be  inserted  in  the  current  sequence, 
if  such  a  node  exists.  Go  to  step  7. 
(The  criteria  by  which  we  determine  if  an 
insertion  can  be  made  are  given  below.) 

Step  5  -  If  there  is  no  such  place  to  insert  node  k, 
then  try  to  find  a  node  L  (ISWAP)  in  the 
current  sequence  such  that  k  can  replace 
I (ISWAP)  and  I  (ISWAP)  has  a  good  cnance  of 
being  reinserted  somewnere  else. 
Select  ISWAP  which  has  the  largest  time 
window  width  among  candidates  for  ISWAP. 
If  there  is  no  candidate,  then  stop. 
<  ■  no  feasible  solution  ■  ) 

Step  6  -  Do  swap  (  add  k  to  N»,  and  delete  L(ISWAP) 
from  N1,  and  set  k  =  L(ISWAP)  ), 
then  update  slack  and  arrival  times. 
Go  to  step  3. 
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Step    7    .    Select    the   node  which   results    in    the    minimum 
additional   travel,      i    .  e,      the      node   k    which 
minimizes    the    following    quantity. 
c(L(I),k)     +c<k,L(I+1))    -   c(L(I)  ,L(1+1))  . 

Step   8    .    Insert    k  after    L(ISTAfi)  ,    and   add    k      to   N', 
and  update   slack  and    arrival  times. 

Step    9    .    If   N»=N,    then    stop. 

Otherwise,    go    to   step    3. 

End   of   algorithm   SLACK 

This       procedure  starts         with     sorting      an      array 

u      ,    u    ,...,u        into        ascending      order      using      a      heapsort 
1  2  n 

[Ref.       20].  This   u    array      is    used    to      select    a   node      k   in 

ascending  order  for  insertion.  Since  the  upper  bound  cannot 
be  violated,  this  step  is  performed.  Then  find  a  node  1{I) 
after  which  ncde  k  can  be  inserted  in  the  current  sequence, 
if    a  such   a    ncde   exists. 

There   are      two   tests    which      must    be      administered   to 
determine  if      k   can      be  inserted      after    L(I).  First,      the 

arrival  time  at  node  k  if  k  succeeds  L(I),  whicn  is  called 
TEST1  must  not  be  greater  than  the  upper  bound  u  .  Second, 
if  k  precedes  L  (1+  1)  ,  tnen  the  resulting  delay  in  arrival  at 
1(1*1),  which  is  called  TEST2  ,  must  not  greater  than 
SIACK(L  (1+1)  )  .      We   can   calculate    TESI1,    TEST2    as   follows. 

TEST1    =  Arrival    time   at    node  k    if    k    succeed   L(i). 

=  max    {AERVI(L  (i)  ),    l(I(i))}    +    d  (L  (i)  )     *■    c(L(i),k). 

TEST2   =  Delay    in    arrival   at    L(i+1)    if   k   precedes   L(i+1). 

=  max    (    TEST1,    1  (k)     }    +    d  (k  )     +    c  (k,  L  (i+ 1)  )  . 
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If  there  exists  more  than  one  node  L(I)   after  which 

node  k   can  be  inserted,   we   select  L  (I)  according   to  the 

criterion  of  least  additional  travel  time-  This  additional 
travel  time,  called  TEST3,  is  given  by 

TEST3  =  C(IjI),k)  +  c(k, 1(1  +  1))  +  C(1(I)  ,L(H-1)  )  . 

When  we  insert  node  k  after  L(I)  ,  we  update  the 
arrival   times      and   SLACKS.  In    the      updating    process,         we 

compute  updated  values  of  SLACK  only  for  the  nodes  whose 
SLACK   actually    changes   as   a    result   of    the   insertion. 

If  there  is  no  place  to  insert  node  k,  we  call  a 
sunroutine  called  'TSWAP'.  TSWAP  tries  to  find  a  node 
L  (ISWAP)  in  the  current  seguence  such  that  k  can  replace 
L  (ISWAP)  and  L(ISaAP)  has  a  good  chance  of  being  reinserted 
somewhere  else.  TSWAP  uses  TEST1  ans  1EST2  to  find  a  candi- 
date for  ISWAE  and  then  uses  a  largest  time  window  width  to 
select  ISWAP.  If  there  exists  such  a  ISWAP,  then  we  do  the 
swap  and  update  SLACKS  and  arrival  times  and  try  to  insert 
again. 

C-       EXACT    SOLUTION    TECHIIQOES    FOB    HABD    TIME    WINDOHS 
1.       Stare-Space   Relaxation   Procedure 

A  dynamic  programming  model  of  the  time-constrained 
TSP    has    been    developed    by    Christof ides    et    al.    [Ref.    5].  We 

applied  their  approach  to  compute  bounding  information 
within    a  branch    and   bound   algorithm. 

Consider  the  TSP  defined  on  the  graph  G  =  {li,\}  with 
the  time  window  constraints,  where  N  is  a  set  of  all  nodes 
of  G,  and  A  is  a  set  of  arcs.  Let  R{j)  be  the  set  of  nodes 
from  which  it  is  possible  to  go  directly  to  node  j.  We  can 
initially  set  R(j)  =  N  -  [  i  |  1;  +  d-L  +  cr  >  u j  } ,  because 
it  is  impossible  to  go  directly  from  ncde  i  to  node  j  if  the 
earlist  possible  arrival  time  at  node  j  exceeds  the  upper 
time   window    cf   node    j. 
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Let  f <S,j)  be  the  duration  of  the  least  time  path 
starting  at  node  1  passing  through  every  node  of  SSS'=N-{1} 
and  finishing  at  node  j.  For  a  given  S  and  j ,  we  can  calcu- 
late a    minimua   arrival   time    in   node    j    as 

T(S,j)    =.      min     f    f  (S-j,i)    ♦    d.+   c.  .].  (3.1) 


Then, 


f  (S,j)     =       T(S,j)  ,    if       1  .    <    T(S,j)       <    U. 
=      1.  ,    if       TD<S,j)        <    1. 

oo  ,    if       T  (S,j)        >u 


3 


with   the  initialization: 


f  <{j}  ,  j)    =  c  (1,  j)         ,    if      1.    <      c        <    u 
=1  ,    if      c3      <    l'3  3 

=         oo  ,      if        C  >     U 

1j  j 

In  equation  (3.1)  the  minimum  arrival  time  in  node  j  passing 
through  the  nodes  in  the  set  S  can  be  described  as  the  sum 
of  three  terms  :  the  first  is  the  duration  of  the  least  time 
path  passing  through  the  nodes  in  the  set  S- {j}  and  ending 
in  node  i,  the  second  is  the  time  required  to  spend  in  node 
i,    and    the    third    is   the   travel   time    from   node    i    to   node    j- 

The  f  (S,j)  can  be  calculated  for  all  subsets  3  of  S« 
and  for  all  nodes  j  by  using  equation  (3.1)  recursively. 
Finally,    the    optimum   solution   can   be   calculated    as 

min      [     f<S,,i)     +    d      +    c         ]. 
i£S'  i  i1 

Since   the   computer  storage   requirements   increase 

exponentially  with  the  size  of   the  problem,   tnis  method  is 

limited  to  small  problems.   The  total  number  of  f(S,j),  when 

'n-l 


S   contains    k    nodes,      is  k(    ,,     I,      since    f(S,j)       must    be   calcu- 
lated  for   all   subsets   S   oi   S',    and    since   each    node    in    S    oust 
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be  considered  as  a  possible  end-node  j.  Therefore  the 
storage  requirement  for  f  (S, j)  in  a  n  node  problem,  is  given 
by  [Eef .    21  ]. 


n-1         /n-1\ 


n-2 
=    (n-1)     2    .  (3,2) 


The  storage  requirements  to  solve  a  22  node  problem 
exceed        22,020,096.  For        relaxing        this         limitation, 

Christofides  et  al.  [fief-  5]  proposed  a  state  space  relax- 
ation procedure  which  is  analogous  to  Lagrangean  relaxation 
[Ref.  22]  in  integer  programming.  The  state  space  associ- 
ated with  a  given  dynamic  programming  recursion  is  relaxed 
in  such  a  way  that  the  solution  to  the  relaxed  recursion 
provides  a  bound  which  could  be  embedded  in  a  general  branch 
and  bound  method  [  Ref .  23].  We  describe  Christofides  et  ai's 
method    for    doing   this    below. 

Consider  the  dynamic  programming  formulation  (3.1) 
The  state  variable  in  that  formulation  is  (S, j) ,  and  the 
stage  is  the  cardinality  of  S.  Let  g  (S)  be  a  mapping  from 
the  domain  of  (S,j)  to  some  other  vector  space  (g(S),j)* 
let  : 

H(g(S),j)    =    {    (g(s-j),i)    |    i  e  (s-jn  R(j))    }  (3.3) 

Since  we  are  interested  in  lower  bounds  to  the  TSP 
with  time  constraints,  H(g(S),j)  in  (3.3)  may  be  replaced  by 
any  larger  set  that  is  easier  to  compute.  Thus,  H(g(S),j) 
can    be    defined    by   the    following   equation: 

H(g(S),j)     =     [    (g(S-j),i)     |     i£E(g(S),j)     }  (3.4) 

where       (S-  j  D  E  (  j) )  £  E  (g  (S  )  ,  j  )  . 

For    calculating      the    lower  bound      of, the      problem,      equation 
(3.1)    can   be    changed    to    the    following    equation: 

T(q(S),j)    =         min        ,  #[f  (g(S-j)  ,i)    +   d.    +  c.  .]      (3.5) 
(g(S-:)  ,i)  e  H  (g(ST,:j)  i  JO 
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min  [f  (g(S-j)  ,i)    ♦    d      +    c      ]       (3-6) 

i  e  E  (g(S),  j)  i  ij 


This   gives    us: 


f<g(S),j)       =      T(g(S),j),   if   1  .<    T(g(S),j)     <    u 


=      1 


,   if  T?g(S),j)    <   1  . 


,   if  T  (g(S)  ,j)    >   u    . 

3 


With   the   initialization; 


f  (g(j)  ,  j)    =     c  ,      if      1.    <    c        <   u. 

1D  3         JD,         D 


=      1  ,      if     c  <      i 

j  ..        Ij  .         J 

,      if     c  >      u 

1j         j 


oo 


Finally,    tne    optimum   solution   can   be   calculated    as 

min  [    f  (g  (S1  )  ri)    +   d        +    c         ]- 

iEE(g(S),1)  i  i1 

The  mapping  can  be  selected  frcm  any  separable  func- 
tion. Christofides  et  al.  used  the  following  mapping 
function. 

g  (S)        =    |S|.  (3.7) 

Then  eguation    (3.6)    becomes    : 

T(k,j)     =      min  [     f(    k-1,i    )     ♦    d      +    c         ]  (3.8) 

i  e  E(k,  j)  i  ij 

where      k   =    | S|    >    1. 

This   gives    us: 

f  (k,j)     =    l(k,j)  ,   if   1<    T(k,j)    <    u 

=   1  .  ,   if   Tik,j)    <   1  .        3 

=    oo  ,   if   1  (k, j)    >   u   . 

3 

With   the   initialization: 

f  dr  j)    =      c(1,  j)       ,      if      1      <   c        <   u 

=  1       ,  if  Cj  <  1i    j 

j  .  ,         1j    .  3 

=         oo  ,         1  f         C  >         U 
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Finally,   the    optimum   solution   can   be   calculated    as 

min  [    f  (JS» |  ,i)     +   d         ♦   c         ]. 

i  £  E(  |N|,  1)  i  i1 


2.      Additional   Condition 

In  the  previous  section,  we  discussed  Christofides 
et  al. • s  state  space  relaxation  procedure  which  provides  a 
lcwer  bound  on  the  TSP  by  reducing  a  state  space  in  dynamic 
programming-  This  lower  bound  is  effective  in  branch  and 
bound  only  if  it  is  a  tight  hound.  This  is  similar  to  the 
case  in  integer  programming  where  the  effectiveness  of 
Lagrangean  relaxation  in  producing  bounds  is  relative  to  tne 
integer  programming  formulation.  A  redundant  state-space 
condition  can  be  helpful  to  get  a  better  bound-  For  this 
purpose,  an  additional  condition  was  used  by  Christofides  et 
al.  to  avoid  loops  formed  by  three  consecutive  nodes  [Ref. 
5].      This  can   be   done    in    the    following    way. 

let  k  =  | S 1 .  let  f  (k,j,1)  be  the  duration  of  the 
least  time  path  from  tne  initial  state  to  state  (k,j) 
without      loops      formed      by   three      consecutive      nodes.  let 

f  (krj,2)  be  the  duration  of  the  second  least  time  path  from 
the  initial  state  to  state  (krj)  without  loops  formed  by 
three  consecutive  nodes.  let  p(k#j,m)  be  the  predecessor  of 
j  on  the  path  corresponding  to  f(k,j,m).  With  the  above 
definition,    recursion    (3.8)     becomes: 

T(k,j,  1)    =.    min       .[     f(x-1,i,m)     ♦    d.+   c.       J,  (3.9) 

i  e  E(k,]J  i        13 

where  m  =  1,  if  p(k-1#i, 1)*j 
=  2,  otherwise. 

This  gives  us: 

f(k,j,1)  =  I(k,j,1),   if  1.  <  T(k,j,1)  <  u.       (3.10) 

=  1.      ,   if  Tik,j,1)  <  1. 
3  3 
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=     00  ,       if    T(k,j,1)     >    u    . 

J 

Recursion     for     f(k,  j,2)         can      be      written      in      the 
following   way. 


Let 

T(k»j,2)    =    min  [    f(k-1,i,m)     +    c 

1  ID 


T(k,j,2)    =    min  [    f(k-1,i,m)     +    d     +    c        ],  (3.11) 

.  i  e  E  <Jc  ,  jl 
i   *p(k,3,1) 

where   m   =    1,    if    p(k-1,i,1)*j 

=   2,    otherwise. 

This  gives    us: 

f(k,j,2)    =    T(k,j,2),         if      1.    <   T(k,j,2)    <    u.  (3.12) 

=   1.  ,        if     T^k#j,2)    <   1. 

=    oo  ,         if      T(k,j,2)    >    u   . 

J 

The    initialization    is 

f(1,i,1)    =    c(1,i)       ,    if      1      <      c        <    u  (3.13) 

i  1i  i 

=      1  ,    if      c  <   1 

i  1i  i 

,    if      c  >    u 

1i  i 


=  00 


and 

f(1,i,2)    =      =o  (3.14) 

Finally,    the    optimum   solution   can  be   calculated    as 

Bin         [    £(|S'|fifD    ♦    d      +    c         ].  (3.15) 

i  e  E(J  N|  ,1  )  i  i1 


Since  the  additional  condition  can  avoid  considera- 
tion of  a  useful  lower  boundr  we  considered  f  (k-1#i,2)  in 
recursion  (3.9)  and  (3.11)  only  when  the  predecessor  of.  i  on 
the  path  corresponding  to  f(k~1*i*1)  is  j.  If  we  do  not 
consider  the  second  least  time  path  in  case  of  p  (k-  1,  i,  1)  =  j, 
then    f(g(S)rj)    dees   not  guarantee    the    lower  bound    of    f(S,j). 

For  this  example,  let's  consider  a  4  node  TSP  with 
time  constraints.  Node  A  is  the  starting  node.  D  is  the 
time  free  node.  The  lower  bound  of  node  B  is  9,  the  upper 
bound   of  node   B    is    11,    tne   lower    bound   of   node    C    is    19,       and 
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tie  upper  boucd  of  node  C  is  21.  Suppose  service  time  at 
each  node  is  zero.  Figure  3.7  shows  an  optimal  route  for 
this    problem. 


Figure    3.7        Optimal   Route   of    Four   Nodes   Problem. 


From   equation    (3.13)    we   can    get: 
f  (1,B,  1)     =    10, 

f  (1,C,1)    =   19, 

f  (1f  Df  1)    =    7.07 

Now    applying    equation    (3.9)     recursively   with   i=1,    for   k  =  2    we 
can    get: 


52 


f  <2,B,1)   =   oo  r 

f <2,C, 1)  =  19, 
f  (2,D,1)  =  17.07 

Similarly,  for  k=3 

f  <3,B,  1)  =  -  , 

f  (3,C,1)  =~  , 

f  (3,D,1)  =~  - 

We  can   see  easily  that   f(3#D,1)   is   not  a  lower   bound  of 
f  (iB.C.L},Q). 

3-   Branch  and  Bound  Procedure 

In  this  section  we  introduce  brancn  and  bound 
enumeration  which  is  used  to  eliminate  subtours  in  the  solu- 
tion of  the  state  space  relaxation  procedure-  Since  the 
state  space  relaxation  procedure  is  a  relaxation  of  the  TSP 
with  time  constraints,  the  solution  to  tne  state  space 
relaxation  procedure  provides  a  lower  bound  on  the  optimal 
value  of  the  ISP  with  time  constraints.  Any  heuristic  solu- 
tion can  provide  an  upper  bound.  He  denote  some  notation  to 
explain  this  algorithm  as  follows. 

FL3D  =  The  lower  bound,  which  is  tne  optimal  solut- 
ion to  the  state  space  relaxation  procedure, 
on  the  optimal  solution  to  the  TSP  with 
time  constraints  given  restrictions  at  the 
current  node. 

Z     =  Current  upper  bound. 

STACK  =  Array  which  represent  decision  tree.  It  con- 
tains arc  lists  which  have  the  same  head  in 
optimal  tour  to  the  state  space  relaxation 
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procedure  given  restrictions   at   the   current 

node. 

£c»  ]  =  Travel  time  matrix  given  restrictions  at 

ij 

the   current    node. 

There  are  two  types  cf  tree  search.  One  is  depth- 
first  search,  the  other  is  breadth-first  search  [Ref.  24]. 
Ke  used  depth-first  search  since  breadth- first  search 
reguired  substantially  more  storage.  Depth-first  search 
simply  means  that  when  a  separation  is  defined,  one  of  the 
nodes  created  by  the  separation  is  immediately  selected  to 
be  the  next  subproblem,  and  when  a  ncde  is  fathomed,  the 
enumeration  always  backtracks  to  the  most  recently  created 
live   node. 

One  of  the  most  important  reguirements  of  any  branch 
and  bound  algorithm  is  tight  bounds.  The  closer  the  bounds 
are  to  the  optimal  solution,  the  fewer  nodes  must  be  enumer- 
ated. We  used  the  SCCO  heuristic,  which  was  described  in 
section  B.2,  as  an  initial  upper  bound.  The  lower  bound  is 
obtained  from   eguation    (3.15). 

To  save  computing  time  we  need  a  criterion  to  decide 
whether  or  not  the  branching  should  be  continued.  If  FLBD 
is  greater  than  Z,  then  the  node  is  fathomed  since  explicit 
enumeration  need  not  be  extended  belcw  the  current  node. 
For  branching  we  consider  the  arcs  which  have  the  same  head 
node  in  the  directed  graph  since  each  arc  must  have  a 
different    head      in    the   TSP    solution.  If    there    is      no    such 

arc,  then  that  solution  is  a  feasible  solution.  After  all 
nodes  of  the  tree  are  fathomed,  a  feasible  solution  which 
has  the  same  value  as  the  upper  bound  is  an  exact  solution 
to    the    TSP    with    time    window    constraints. 

The  following  branch  and  bound  algorithm  is  used  in 
the    programs    written   for    exact   solution- 


54 


Algorithm  :  Eranch  and  Bound  Procedure 

Input  :  Total  travel  time  of  heuristic,  travel  time. 
Cutput:  Ordered  list  of  tour,  total  travel  time. 

Step  1  .  (Initialization) 

Let  Z  =  the  optimal  solution  of  SCCO. 

STACK  =  empty. 

[    c'.]  =   [    c.       ] 
Step    2    .    Compute   FLBD    given    restrictions   defined    by 


[    c»     ].    If   FLED  >    Z,    go    to    step    5. 
ij 
£tep    3    .     (Construct    the  tree) 

Put    all   arc(i,j)     which    have    the    same   head    j 

in   directed    graph   on    STACK. 

If   there    is    no   such    arc,    save    feasible    route 

and   update   Z    =   FL3D    then   go    to    step   5- 
Step    4    .    Let    travel  time   of    arc(i,j)    which   is   in    the 

top   of   STACK    be   infinite,    then    go    to   step    2. 

(i.e.  ,    c  '      =  °°  . ) 
Step    5    .     (Backtrack) 

If    STACK    =   empty,    go    tc    step    7. 
Step    6    .    If    travel   time  of   arc(i,j)    which    is   in    the 

top  of   STACK    is    finite,    let    travel   time    of 

that      arc  (1, j)      be      infinite,       then   go    to 

step    2.     (i.e.,    c*         =  °°   -) 

ij 
Otherwise,    let  travel   time   of    that    arc(i,j) 

be   original   travel    time    of    that   arc(i,j)    and 

remove    tnat    arc(i,j)    from   top   of    the   STACK. 

Go    to   step  5.     (i.e.  ,    c*      =    c      ) 

ij  i  j 

Step    7    .     (termination) 

If    there      is    a      feasible   route,    then   the 

optimal    travel   time    =    Z. 

Otherwise,    there    is    no    feasible    solution. 

End  of   algorithm   Branch    and    Bound    Procedure 
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We    present   the   results      of    our   computational  experi- 
ence   with   the   algorithms    of    this   Chapter   in    Chapter   V. 
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IV.  THE  ISP  WITH  SOFT  TIME  ilJLQCi  CONSTRAINTS 

A.   INTBODUCTION 

The  second  time-constrain ted  TSP  we  consider  is  the  case 
in  which  both  late  and  early  arrivals  are  allowed  by  paying 
a  penalty  cost.  The  penalties  are  allowed  to  be  different 
for  early  and  late  arrivals.  Tne  penalty  cost  is  calculated 
as  follows. 

Upper  penalty  cost  =   max  [  0,   upper  penalty  constant 

x  {  arrival  time  -  upper  bound  )  ]- 

Lower  penalty  cost  =   max  [  0,   lower  penalty  constant 

x  (  lower  bound  -  arrival  time  )  ]- 

In  fact,  the  upper  penalty  constant  is  greater  than  the 
lower  penalty  constant  in  most  cases.  Figure  4.1  may  be 
helpful  to  understand  this  case. 

This  approach  makes  every  problem  feasible,  no  matter 
what  the  time  windows  are,  i.e,  even  if  it  is  infeasible  in 
the  hard  time  window  case.  This  reflects  a  practical  point 
of  view,  especially  when  it  is  possible  to  save  a  great  deal 
of  mileage  by  allowing  a  small  amount  of  time  window 
violation. 

In  this  Chapter,    we  considered  one  unit  of   cost  to  be 

the  same  as  one  unit  of  time.   In  real  world  problems,  it  is 

possible  to  get  a  cost  by  multiplying  traveling  time  by  some 

constant. 

We  use  the  notation   lp   and   up   for   the   lower   and 

k'        k 
upper  penalty  cost  at  node  k. 
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early 
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Figure    4. 1        Diagram   for   Soft    Tine   Window    Case. 

B-       HEURISTIC    S010TIOB   TECHSIQUES    FOB    SOFT    TIME    IINDOHS 

1-       Nearest    Neighbor 

This  heuristic  is  similar  to  the  hard  time  windows 
except  it  takes  into  account  any  penalty  cost  that  might  be 
necessary. 

Algorithm    :    nearest   Neighbor 

Input       :    Number    of    nodes,    x  and    y   co-ordinates  of    all 

nodes,    time   windows    for    all    nodes. 
Output    :   Ordered   list    cf    tour,    tctal   cost. 

Step    1     .    (Initialization) 

Start    at    the    depot. 

Let    i=depotr    N*    =    {i} ,    cost    =   0- 

i 

Step    2    .    Compute   AftRVT    tor   all-ncdes    k€N-N' 

AEEVT    =    AEEVT       +    d       ♦    C       . 
k  i  i  ik 

cost      =   cost       +   d      +    c 
k  i  i  ik 


58 


Step    3    -    If   ARRVT      <   1       ,    then   cost   =   cost    +    lp      . 

Jc  Jc  k  k  k 

If   AERVT      >   u       ,    then   cost   =  cost   «•    up 
k  k  k  k  k 

Step    4    .    (  Nearest    Neighbor   Selection  ) 

Select    the    node   k    €    N-N •    such   that   cost 

k 
is  a   minimum.     I.e.    ,    find 

k  =    argmin      cost  s.t.    j  6N-N". 

J 

Step  5  .  (  Insertion  ) 

Insert  k   after  i,  add  k  to  subtour  N',  and 

let  i  =  k  . 
Step  6  .  If  N'  =  N  ,  then  go  to  next  step. 

Otherwise,  go  tc  step  2. 

Step  7  .  Compute  total  cost,  then  stop. 

Total  cost  =  ccst  +  d  +  c 

k    k   k, depot 

End  of  algorithm  nearest  Neighbor 

This  solution  was  constructed  by  starting  at  the 
depot  and  moving  to  the  nearest  neighboring  customer  that 
has  not  yet  been  visited.  The  term  "nearest"  is  modified  in 
the  sense  tnat  we  add  a  penalty  cost  to  the  travel  time  if 
the  time  window  for  city  i  is  violated. 

2.   SCCO 

This  algorithm  is  also  designed  for  the  case  when 
there  is  a  combination  of  tight  time  window  nodes  and  time 
free  nodes.  The  strict  observance  of  the  upper  bound  in  the 
hard  time  windows   is   replaced   by   a    penalty   cost. 

Algorithm    :    SCCO 

Input       :    Number    of    nodes,    x   and    y    co-ordinates  of    aj.1 

nodes,    time   windows    for   all    nodes. 
Output    :    Ordered    list    of   tour,    total    cost. 

Step    1    .     (Initialization) 
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Start  at  the  depot. 

let  i=depot,  N'  =  {i}  ,  cost  =  0. 

i 

Step  2.   Set  k  =   argmin   u    s.t.   j€N-N*. 

If  k  is  time  free  node,  then  set  k  =  depot. 

Step  3  -  Insert  node  k  in  the  subtour  N'. 

Compute  ARRVT 

k 

ARRVT   =  ARRVT   ♦  d    +  C 

k        i    i      ik 

Step  4  .  Insert  time  free  node  j€  N-N '   between  nodes 

i  and  k   by   cheapest  insertion  and  cheapest 

selection  (  same  as  CCCC)  until  ARRVT    does 

k 

not  exceed   u  . 
k 

Step    5    .    Update    cost 

k 

cost   =      cost    ♦      d      +    c 
k  i  i  ik 

If  ARRVT  <  1  ,  then  cost  =  cost  +  lp 
k    k  k      k     k 

If  ARRVT  >  u  ,  then  cost  =  cost  +  up 
k   k  k      k    "  k 


Step  6  -  let  i  =  k. 

If  N'  =  N,  then  go  to  next  step. 
Otherwise,  go  to  step  2. 

Step  7  .  Apply  tne  Modified  Oropt  procedure  to  the 

current  tour.  Stop  when  no  further  improve- 
ments can  be  found. 

End  of  algorithm  SCCO 

This  procedure  is  also  similar  to  the  cheapest 
selection,  cheapest  insertion  method  fcr  the  unconstrained 
TSP,  except  that  the  nodes  with  time  windows  are  treated 
differently  froc  the  time  free   nodes.    The  nodes  with  time 
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windows  are  inserted  in  order  of  increasing  upper  time 
window  bounds.  The  time  free  nodes  are  inserted  between 
those  nodes  ty  cheapest  selection  and  cheapest  insertion, 
for  as  long  as  the  upper  bound  of  the  time  windows  will 
allow. 

In  the  end,  a  Modified  Oropt  is  used  to  improve  the 
solution.  This  procedure  consider  only  those  exchanges  that 
would  result  in  a  node  being  inserted  netween  two  other 
nodes   in  the    current    tour. 

3.       SCAO 

This  algorithm  is  also  designed  for  the  time  window 
set  which  is  composed  of  some  tight  time  windows  and  some 
time   free  nodes. 

algorithm    ;    SCAO 

Input       :    Number    of    nodes,    x   and    y    co-ordinates   of    all 

nodes,    time   windows   for   all   nodes. 
Output    :    Ordered    list    of   tour,    total    cost. 

Step    1    .     (Initialization) 

Start   at    the    depot. 

Let    i=depot,    N1    =     {i},    cost    =    0. 

i 

Step    2.      Set   k   =      argmin      u  s.t.       j£N-N*. 

3 

If   k   is   time    free    node,    set    k   =   depot. 

Step    3    .    Insert    node  k    in   the    subtour    N'. 

Compute    AfiEVT 
k 

AEEVT      =    ABBVT       ♦    d         +    C 

k  i  i  ik 

Step   4    .    Insert    time  free  node    j    €   N-N*    between   nodes 

i   and   k    by      cheapest    insertion    ana      greatest 

angle       (same    as   CCAO)  until    AfiEVT    does    not 

k 

exceed      u    . 


b1 


Step  5  -  Update  cost 

k 

cost  =   cost  ♦   d   +  c 
k       i     i    ik 

If  ARRVT  <  1   r  then  cost  =  cost  +  lp 
k   k  k      k    k 

If  ARRVT  >  u  .  then  cost  =  cost  +  up 
k   k  k      k     k 

Step  6  -  let  i  =  k. 

If   N'    =   N,    then  go   to    next   stop- 
Otherwise,    go    to  step    2. 

Step    7    .    Apply    the    Modified-Orop t    procedure    to   the 

current  tour.     Stop   when   no   further   improve- 
ments  can    be    found. 

End  of   algorithn   SCAO 

This  algorithm  is  same  as  SCCO  except  a  greatest 
angle   selection    in   stead    of    a    cheapest    selection    in   SCCC. 

C-       EXACT    SOIOTICM    TECHHIQOES    FOB    SOFT    IIME    IINDOIS 
1.       State-Space  Relaxation   Procedure 

In  this  section  we  describe  a  state  space  relaxation 
procedure,  which  is  adapted  from  Christofides  et  al.  [Ref. 
5  ]#  for  soft  time  windows.  They  only  considered  the  TSP  with 
hard  time  windows  and  without  time  windows.  The  difterences 
are  as  follows.  Tne  waiting  cost  is  replaced  by  a  penalty 
cost  to  be  paid  in  the  early  arrival  case.  Late  arrival  is 
allowed,  but  a  penalty  cost  has  to  be  paid.  So  we  have  to 
calculate  tne  duration  and  the  penalty  cost  on  each  possible 
path  to  decide  the  least  cost  path  in  eacn  stage.  Be  denote 
the  penalty  ccst  on  each  possible  path  as  PC  in  this 
section . 

Consider  the  TSP  defined  on  the  graph  G  =  {N,A}  with 
soft  tine  window  constraints.  Let  S*  be  a  set  of  all  nodes 
except    starting    node.       Let  S    be   a   subset   of    S'.         Let   f  (S, j) 
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be  the  cost  of  the  least  cost  path  starting  at  node  1 
passing  through  every  node  of  S  and  finishing  at  node  j. 
Let  T(S,j)  be  the  total  duration  of  a  path  corresponding  to 
f(S,j).  Let  p(S#j)  be  the  predecessor  of  j  on  the  path 
corresponding   to  f <S, j) -  Let   lp(t)      be  the      early  arrival 

penalty  cost  function  and  up (t)  be  the  late  arrival  penalty 
cost  function-  For  a  given  S  and  j,  total  duration  of  a 
path  can  be    calculated   as 

T  (S,j)    =   [    T(S-j,i)    +    d.    ♦   c.       ].  (4.1) 

where   p  (S,  j)     =    i- 

Id  eguation  (4.1)  total  duration  of  the  least  cost  path 
passing  through  the  nodes  in  the  set  S  and  ending  in  node  j 
can  be  described  as  tne  sum  of  three  terms:  the  first  is 
total  duration  of  the  least  cost  path  passing  through  the 
nodes  in  the  set  S-{j}  and  ending  in  node  i,  the  second  is 
the  time  required  to  spend  in  node  i,  and  the  third  is  the 
travel  time  from  node  i  to  node  j.  The  dynamic  programming 
recursion  to  determine  the  least  cost  path  may  then  be 
stated    as 

f(S,j)     =    min      [    f(S-j,i)    ♦    d    ♦    c         *    PC    j  (4-2) 

i  £  S-  j  i         i  j 

where  T1   =    [    2(S-;j,i)    *    d      ♦    c         ]. 

i  ij 

PC=0  ,if      1      <  T  1      <    u 

=    lp  (1    -T1)     ,    if      T  1       <    I3 


=    up  (Tl-u    )     ,    if      11       >    u 


with   the  initialization: 

f  ({j}  r  j)     =    C  ,     if        1        <        C  <     U 

1J  j  1 j  j 

=   c         ♦   lp  (1    -c      )     ,     if      c         <    1 
=  .  c        +    up  (c      -u    )    ,    if      c        >    u 

Finally,  the  optimum  solution  can  be  calculated  as 

min   [  f(S,,i)  *  d   +  c    ]. 
i  e  S»  l      i1 
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Since  the  computer  storage  requirements  are 
increased  exponentially  with  the  size  of  the  problem,  this 
method    is      limited    to      small    problems.  For   relaxing      this 

limitation,  a  state  space  relaxation  procedure  can  be  used 
same  as    Chapter    III. 

Consider  the  dynamic  programming  formulation  <4.2) 
The  state  variable  in  that  formulation  is  (S,j),  and  the 
stage  is  the  cardinality  of  S.  Let  g  (S)  be  a  mapping  from 
the  domain  cf  (S,j)  to  some  other  vector  space  (g(S),j). 
Let: 

H(g(S),j)    =    {    (g(S-j).i)    |    i€S-j}  (4.3) 

Since  we  are  interested  in  lower  bounds  to  the  TSP 
with  time  constraints,  H(g(S),j)  in  (4-3)  may  be  replaced  by 
any  larger  set  that  is  easier  to  compute.  Thus,  H(g(S),j) 
can    be      defined    by    the   following    equation; 

H(g(S),j)     =    {    (g(S-j),i)     I    i€E(g(S),j)     }  (4. 4) 

where       £-j  £  E  (g  (S)  ,  j)  . 

For  calculatirg  the  lower  bound  of  the  problem,  recursion 
(4.1)    can    be    changed    to   the    following    equation; 

T  (g(S)  ,  j)     =   [     T(g(S-j)  ,i>     ♦    d.     +    c.        ]  (4.5) 

where      p(g(S)  ,j)    =   i. 

Recursion    (4.2)     may    be    stated    as 

f(g(S),j)    =.      min        [f  (g(S-j)  ,i)+    d.+    c.       ♦    PC]       (4.6) 
(g(S-^)  ,i)£    H(g(Sf#3)  r        lj 

min        £f  (g(S-j)  ,i)  «■    d   +   c        +    PC]       (4.7) 
i  e  E  (g(S),  j)  i        ij 

where    T 1   =   [     1(g(S-j)  ,i)     +   d    +   c         ], 

PC   =    0  ,if      1      <   T  1      <    u  t    - 

=    lp  (1  .-T1)     ,    if      T  1       <    I3 

=    up  (Tl-u    )     ,    if      T  1       >    u 
j  j 


64 


with  the  initialization: 

f  (g(j)  ,  J)    =   c    .  ,    if      1.    <  c  <   u. 

=  c          ♦  lp(l  -c      j3      ,  if3  c      3<    1 

1j  #    J     1j    t  .  _  1]s      J 

=  c          *  up(c        -a    )     ,  if  c        >    u 

1j  1j        J  1j           J 


Finally,    the    optimum    solution   can   be   calculated    as 

min  [    f  (g(S*  )#i)    ♦  d        ♦    c        ]. 

ie£(q(S),1i  i  i1 


The  mapping  can  be  selected  frcm  any  separable  func- 
tion. He  used  a  mapping  function  (3.7)  ,  which  is  proposed  by 
Christofides    et    al. ,         same   as  Chapter    III.  Then    equation 

(4.  5)    becomes: 

T  (|S  |  ,  j)    =   [     T(|S|-1,i)     +    d.     +    c.  .     ].  (4.8) 

i  lj 

where      p(|S  |,  j)     =   i 

Recursion    (4.7)     may   be  stated   as: 

f  (|S|,j)    =      min      [  f  (|S  |-1,i)+    d      +    c        +    PC  ]  ('4.9) 

i£E(JSj,j)  i  ij 

where   T1   =    [     I(JS|-1,i)    +    d      +   c         ], 

i  ij 

PC  =   0  r  if      1      <  T  1      <u 

=    lp  (1    -T1)     ,   If      T1      <    i3 
=    up  (TT-u    )     ,    if      T  1      >    u 
with   the  initialization: 

f  0» j)    =  c  ,   if     1      <     c  <   u 

=    c  +    lp  (1  -c      )     ,    if      c         <    1 

1j         '     j     11       _      1j  v    j 

=   c  *    up (c      -u   )     ,    if      c        >    u 

1j  1j     j  1j        j 


Finally,    the    optimum  solution   can  be   calculated    as 

air  [    f  (|S»  |  ,i)     ♦  d        +  _c         ]. 

i  e  E(|N|,1)  i  i1 
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2.   Additional  Condition 

In  th€  previous  section,  we  discussed  a  state  space 
relaxation  procedure  which  is  adapted  from  Christofides  et 
al,[Ref.  5].  That  procedure  provides  a  lower  bound  on  the 
TSP  with  soft  tine  window  constraints.  The  additional 
condition  to  avoid  loops  formed  £y  three  consecutive  nodes 
was  used  to  get  a  better  bound  [ Eef .  5].  This  can  be  done 
in  the  following  way. 

Let  k  =  |S|.  Let  f  (k,j,1)  be  the  cost  of  the  least 
cost  path  from  the  initial  state  to  state  (kr j)  without 
loops  formed  by  three  consecutive  nodes.  Let  f  (k,j,2)  be 
the  cost  of  the  second  least  cost  path  from  the  initial 
state  to  state  (k, j)  without  loops  formed  ny  three  consecu- 
tive nodes.  Let  p(k,j,m)  be  the  predecessor  of  j  on  the 
path  corresponding  to  f{k,j,m).  With  the  above  definition, 
eguation  (4.8)  becomes: 

1<krj,m')  =  £  T(k-1,i,m)  +  d  +  c    ],  m'=1,2     (<*.10) 

i    ij 

where   p  (k, j , a' )  =  i 

m  =  1  t    if  p(k-1,i,1)*j 

-    2    ,    otherwise. 

With  the  initialization: 
T  (1*  j,  1)  =  c 

and 

T(1,j,2)  =   *  . 

Recursion  for  f  (k,j,1)    can  be  calculated  in   the  following 
way.   Let: 

T'  (k#j,m)  =  £T(K-1,i,m)  +  d  ♦  c    ],  m  =  1,2 

i    ij 

This  gives  us: 
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f(k,j,1)    =      min        rf(k-1,i,m)     ♦    d     ♦   c        +    PC  ]      (4.11) 
is  E(k,jJ  i  13 

where   PC  =    0  ,    if      1.    <  T»  (k,j,m)       <    u. 

=    lp<l   -TMk.J»«)).    ^      TM^^m?      <    1. 

=   up  <T^(k,j,m)-u  ),    if     T«(k,  j,m)      >   u 

3  J 

in  =    1    ,   if   p  (k-1,i,1)  *j 
=    2    ,    otherwise. 

With   the  initialization: 

f(1,i,1)    =c  ,    if      1      <      c        <   u  (4-12) 

1i  i  1i  1 

=  c    .       ♦  lp(l.-c      )     ,    if     c  <    1. 

1i  1      1i  1i  1 

=   c  +    up(c      -u    )     ,    if      c  >    u 

1i  1i      i  1i  i 

Recursion   for  f  (k,j,2)    can  be    written    in   the   following    way: 

f(k,j,2)    =      min        f  f  (k-1,i,m)     ♦    d      +    c        +    PC]      (4.13) 
i  eEtt.jt  i  ij 

i#P(k,3,1) 

where   PC   =    0  ,   if      1      <  T»  (k,j,m)       <    u 

=   lp(l   -T«  (k,  j?m)),    if      T«(k,  3,111?      <   1. 

=    up|TMk,i,m)-u   ),    if      T«(k,  j,m)      >    u. 

3  3 

0  =    1    r   if  p  (k-1,i ,1)  #j 

=    2    ,    otherwise- 

With   the   initialization: 

f(1,i,2)    =     00  (4.14) 

Finally,    the    optimum   solution    can    be   calculated    as 

fflin         [f  (|S'lriJ)    ♦    d      +    c         ].  (4.15) 

i  e  E  ( |  N  1 , 1 )  i  i  1 

Since  the  additional  condition  can  avoid  considera- 
tion of  a  useful  lower  bound,  we  considered  £  (k-1,i,2)  in 
recursion  (4.11)  and  (4.13)  only  when  the  predecessor  of  i 
on  the  path  corresponding  to  f(k-1ri,1)  is  3.  If  we  do  not 
consider  the  second  least  cost  path  in  case  of  p  (k- 1 ,1,  1) = j , 
then   f (g(S)  , j)    does  not  guarantee   the    lower   bound    of   f  (S,j) 


67 


For  this  example,  let's  consider  4  node  TSP  with 
time  constraints.  Node  A  is  the  starting  node-  D  is  the 
time  free  node.  The  lower  bound  of  node  B  is  9,  the  upper 
tound  of  node  B  is  11,  the  lower  bound  of  node  C  is  19,  and 
the  upper  bound  of  node  C  is  21.  Suppose  service  time  at 
each  node  is  2ero,  lp(t)=t,  and  up(t)  =  5t.  Figure  3.7  shows 
an  optimal  route  for  this  problem.  Frcm  equation  (4.10)  and 
(<J.  12)     we  can   get: 

f  (1,B,1)    =    10,    T(1,B,1)     =    10,     p(1,B,1)     =    A; 

f(1,C,1)    =    19,    T(1,C,1)    =    14.14,     p(1,C,1)     =    A; 

f  (1/D,1)     =    7.07,    T(1,D,1)     =    7.07,    p(1,D,1)     =    A. 

New    applying      equation    (4.10)      and    (4.11)         recursively    with 
i=1,    for  k= 2    we    can   get: 

f  (2,3,1)    =   min    [    94.7,     29-84    ]   =    29.84, 
i  e  {C ,  D} 

T  (2,3,  1)     =    14.14,    p(2,B,1)     =    D. 
Similarly, 

f  (2, C,1)     =    19,    T(2,C,1)    =    14.14,     p(2,C,1)     =    D; 

f(2,D,1)     =    17.07,     1(2, D,1)     =    17.07,    p(2,D,1)     =    3. 

Fcr    k   =   3, 

f  (3,B,1)    =   min    [    94.7     ]   =   94.7, 
i   e{c} 

T(3,B,1)  =   24.14,    p(3,B,1)     =    C. 

Similarly, 

f  (3,C,1)  =    39.84,    1(3, C,1)     =    24.14,    p(3,C,1)     =    D; 

f  (3,0,  1)  =  <*>  . 

We  can      see    easily    tnat      f(3,D,1)       is      not    a    lower      bound   of 
f  ({B,C,C}  ,D)  . 
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3-   Branch  and  Bound  Procedure 

We  used  the  same  branch  and  bound  procedure  used  to 
eliminate  subtours  in  the  solution  of  the  state  space  relax- 
ation procedure  in  Chapter  III.C.  3.  The  SCCO  heuristic, 
which  was  described  in  section  B.2,  was  used  as  an  initial 
upper  bcund,  and  the  lower  bound  was  obtained  from  equation 
(4.15). 

We  present  the  results  of  our  computational  experi- 
ence with  the  algorithms  of  this  Chapter  in  the  next 
Chapter  - 
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V-     COMPUTATIONAL    EXPERIENCE 

A.       TES1    PROBIEHS 

Four  sets  of  test  data  are  used  in  this  thesis-  lest 
problem  number  [1]  is  taken  from  Sedgewick  £Ref.19:  p. 309], 
The  other  problems,  numbered  [2],  [3]  and  [4],  are  from 
Appendix  9.1  cf  Eilon  et  al-  '  s  text  [Bef-  21].  These  test 
problems  are  shown  in  Appendices  A,B,C  respectively.  These 
published  problems  contain  node  and  depot  locations,  but 
they   do   not    include   time    windows. 

We  constructed  time  windows  for  test  problems 
Ml#f2],£3]  by  first  using  the  CCAO  Heuristic  on  the  uncon- 
strained ISP.  lime  windows  were  then  placed  aoout  each  node 
such  that  the  CCAO  route  was  feasible.  The  idea  for  gener- 
ating time  windows  in  this  way  comes  from  Baker  [fief-  25], 
who  used  the  unconstrained  Nearest  Neighbor  heuristic  as  his 
starting   point    instead   of  CCAO. 

The  time  window  widths  were  set  to  varying  sizes  ranging 
from  3  to  14.  Seme  of  tne  time  windows  were  rairly  tight 
while  others  overlapped.  This  is  in  contrast  to  Baker's 
work,    where    all   the   time    windows   have    width   egual    2    units. 

The  last  problem  number  [4]  is  the  same  as  test  prctlea 
[3],  except  that  the  time  windows  were  constructed  from  a 
Nearest  Neighbor  solution  to  the  unconstrained  traveling 
salesman      problem,      as      in   Baker      [  Eef .         25].  Figure      5.1 

displays  the  CCAO  solution  for  test  problem  [3]  and  figure 
5.2  illustrates  the  unconstrained  Nearest  Neighbor  solution 
for  the  test  problem  £4],  We  found  a  small  error  in  Baker's 
TSP  solution  for  the  Nearest  Neighbor  [Ret.  25],  in  that  the 
nearest  node  from  node  16  is  node  17,  not  node  13.  The 
resulting  cost    is   actually   higher,    it    is    312.09,    not    310.22. 
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Figure   5.1        Oncon strained  Solution   Obtained   by  CCAO. 

Each  of  the  four  sets  of  test  data  was  used  to  create 
four  test  problems.  The  separate  instances  differed  in  the 
percentages  of  time  window  constraints  that  were  chosen  to 
be  in  effect-  The  four  cases  were  100%,  90%,  75%,  and  50%. 
We   refer  to    this   percentage   as  the    "time-window    percentage". 

A  random  number  generator  was  used  to  decide  which  nodes 
would  have  time  windows.  Test  problems  for  the  time  window 
constrained    TSJ?   are      shown   in    Appendices   G      through    V.         The 
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>- 


190 


Figure  5.2    Uncon strained  Solution  Obtained  by 
Nearest  Neighbor  Heuristic. 


penalty  cost  factors  can  be  varied  depending  on  real  world 
problems.  He  used  2  and  5  as  the  lower  and  upper  penalty 
cost  factor.  Also  we  set  th a  service  time  at  each  node  to  0 
to  make  it  easy  to  construct  the  time  windows. 

The  computational  results  are  presented  in  Tables  II  and 
III.  The  figures  reported  represent  results  of  our  test 
runs  for  each  test  case. 
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B.   COflFOTATIOHAl  RESULTS 

1-   Hard  lime  Windows 

As  ncted  in  section  3,  the  Nearest  Neighbor 
heuristic  often  cannot  solve  the  problem,  because  the 
arrival  time  of  the  nearest  node  frequently  violates  the 
upper  bound.  However  in  test  problem  [4],  the  results  of 
the  Nearest  Neighbor  are  the  same  as  in  the  unconstrained 
problem,  because  this  problem  itself  was  constructed  by  a 
Nearest  Neighbor  heuristic. 

Generally,  the  SCCO  and  SCAO  heuristic  can  be  easily 
applied  to  the  hard  time  window  TSP.  According  to  our 
experiments,  if  the  time  window  width  becomes  large  relative 
to  the  travel  time  between  nodes  on  the  optimal  uncon- 
strained TSP  route,  then  the  lower  percentage  time  window 
problems  become  more  difficult  to  satisfy.  This  phenomenon 
can  be  seen  in  test  problem  [ 1  ]•  That  is  because  the  other 
nodes  in  the  optimal  route  for  the  unconstrained  TSP  problem 
could  b€  inserted  without  causing  violation  of  the  upper 
bound. 

The  S1ACK  heuristic  takes  slightly  more  time  than 
the  other  heuristics.  It  achieved  lower  accuracy  in  test 
problems  [  1  ]  and  [3]  in  the  50  percent  time  window  case. 
The  exact  algorithm  can  find  the  exact  answer  in  most  prob- 
lems, but  when  there  are  fewer  windows  in  effect,  it  takes 
mere  computation  time.  It  cannot  solve  the  50  percent  time 
window  problems  [2]  and  [4]  within   180  seconds. 
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2.      Soft    lime    Windows 

All  of  the  methods  tested  for  soft  time  windows  were 
able  to  find  some  answer  to  every  problem  within  reasonable 
computing  time,  except  for  two  instances  with  the  exact 
algorithm.  Witn  the  Nearest  Neighbor  heuristic,  the  quality 
of   solution   is  not  desirable.  In    general,      the    lower    time 

window  percentage  problems  have  lower  solution  quality.  As 
in  the  hard  window  case,  on  test  problem  [4],  the  results  of 
the  Nearest  Neighbor  heuristic  coincide  with  the  uncon- 
strained TSP  heuristic,  because  this  problem  itself  was 
constructed    bj   a   Nearest    Neighbor. 

As  in  the  hard  time  window  problem,  SCCO  and  SCAO 
generally  find  an  optimal  solution  except  for  one  problem 
witn  50  percent  time  windows.  In  test  problem  [1]  with  50 
percent  time  windows,  the  SCCO  and  SCAO  values  were  215.686, 
165.544  respectively.  The  exact  algorithm  could  not  solve 
the  two  test  problems  with  50  percent  time  windows  within 
180  seconds.  The  reason  is  that  the  solutions  of  the  state 
space  relaxations  have  many  subtours  and  it  takes  a  long 
time    to   eliminate    these   subtours. 

With  both  hard  and  soft  time  windows,  the  results 
are  sensitive  tc  the  percentage,  width  and  position  of  the 
time   windows.  In   most    problems,         the   fewer      time    windows 

there   are,    the   lower    the    accuracy   of   the   heuristics. 


75 


CO 

EC 
o 

Q 

as 
H 

2* 

M 

so 

H 
EH 

H 
Em 
O 
CO 

H      W 

h    as 
H     Eh 


w 

fc-t 


CM 

O 

CO 

6-1 

B 
CO 

cq 


O 
H 

E-t 

-J 

H 

a 
cu 

* 
O 

u 


a 

I     fno^oo 

coocn 

ooor* 

OOOct 

Oh 

^©^■r-^ 

COr-     t 

C7NOCn\0 

coco  • 

4-» 

CJ 

OOrO 

t-t-^o 

T-fSJ*-       , 

T-r—  CO 

O 

t    *    r   * 

•     ivO 

t    1    If" 

•    tin 

rd 

tt 

^r^^r^ 

cnoo^ 

p-r>>r«>p* 

Cncft'- 

W 

oooo 

CNCNCN  1 

CO  coco  CO 

0OC0CN  1 

-p 

SO  O'nO  O 

ooo 

^■^•=r^ 

OOP0 

CO 

•    •   t    • 

•   i   • 

t  #  •  • 

•    •    • 

o 

vO'-O'-O  o 

CTi<7»CTN 

co  co  oo  ao 

CNCN00 

u 

'^C-O'sO  -0 

'.0-0  gj 

r*r»r-r* 

r-r-O 

• 

^^^r 

CNCNCNCN 

CO  CO  CO 

3 

p»r*-r<»o 

r-*roroo 

f>-OM^ 

ocoor* 

cu 

*-Ot-CN 

Or-*-CN 

0<-r-0 

(Nc\)«N^- 

u 

oooo 

oooo 

oooo 

oooo 

o 

<5 

^^^•^ 

ovo^c^cn 

r^r-r-r- 

<?><J\r-<j\ 

u 

+J 

ooo^ 

CNCNCNCN 

ro  PO^OCO 

COOOCNOO 

CO 

CO 

o 
u 

^Q'-GsOlD 

oooo 

J*^-:*^ 

OOCOO 

^\Q\ais~> 

cncr^OCTi 

CO  co  coco 

CNCN0OCN 

^O^O^O^Q 

^O^OvOvO 

r-r-p-r* 

T-T-O*- 

f» 

^^•^•^ 

(NCNCNCN 

co  ro  CO  00 

r-»r^t^rr> 

r-oooo 

r*»r»or- 

P^OP^po 

o 

00*-r- 

Or\j<— CN 

r-OCsjO 

OCN»-r- 

Cu 

u 

oooo 

OOOO 

oooo 

OOOO 

o 

d-cr=r^o 

CTiCTNCTicn 

r^rr-r* 

O^CTi'- CTi 

u 

QOOCO 

CNCNCNCN 

co  co  co  no 

CO  00  CN  CO 

o 

+J 

^DvflvO^O 

OOOO 

^^■■ir^t 

ooroo 

CO 

co 

•    •    •    • 

till 

•  t   t  i 

#   •   •   • 

o 

vO^O^OLD 

en  en  en  o> 

CO  00  coco 

CNCNC0CN 

u 

^O^O^Or" 

\Ov.0^OvO 

r»r«»r»r^ 

«■■  t—  Or 

CN 

^■^•^■a- 

CNCNCNCN 

CO  co  co  co 

p^r-i^ip* 

oooo 

rOr*»r*co 

COOOO 

a 

0«-*-'~ 

r- CNCNCN 

»-OOCN 

t—  CN«—  r- 

H 

■»->  O 
C0.Q 

CJ 

OOOO 

OOoO 

OOOO 

OOOO 

0JJ3 

vOsOHDo 

crscoov-o 

r^or^^ 

(TiCT>0>C^ 

HOi 

■M 

■     ^ovomLn 

CNCOCNCTt 

cocn^*~ 

00  CO  CO  CO 

(0-H 

0)01  1 
ZZ  1 

o 
u 

r-r^r*«- 

ooo<- 

3-CJNrO^r 

oooo 

C^C^COcn 

OM**  ij*)  CTi 

OOct^^O 

CNCNCNCN 

^&\Qr~  o 

^^3-cncn 

P**-c*o 

«—  *—  r»T— 

OOCN 

=tcno>T~ 

CNvOO^T 

CO  CO  CO  CO 

01 

4J 

0)  3 

u 

a  o  t 

3d> 

m^fNco 

t-o>._o^- 

'—CTiLnr- 

i-omt»«- 

•HT3-r 

H<4-| 

^"r-*- 

(MT-T-T- 

cN»—  *■  «■- 

(N*—  v—i— 

Ci  « 

^ 

•H 

w 

M3 

U>        (. 

'J 

■£»      c 

D 

a        rj 

3 

^^OsCO 

CNCNCNCN 

CNCNCNCN 

CNCNCNCN 

34-1  C 

3 

*— r— t-t— 

CNCNCNCN 

INCNINCN 

CNCNCNCN 

Z  O  c 

3 

a 

CD 

H 

J* 

/-** 

*-»■* 

r^s 

r*« 

o 

r» 

CN 

CO 

^r 

Ul 

*-«< 

>mmi 

<-w 

Wtf 

a, 

CO 

CO 
O 
ro 

CQ 


a, 
o 

CO 

a 
o 
o 

0) 
CO 


CO 

a> 

a 

•H 


o 

Cu 


76 


VI.  COHCLOSIOMS  AMD  RECOMMENDATIONS 

This  thesis  has  presented  some  heuristics  and  exact 
algorithms  for  the  solution  cf  traveling  salesman  problem 
with  tine  window  constraints.  We  considered  two  different 
kinds  of  time  window  constraints  :  hard  time  windows  and 
soft  time  windows.  Hard  time  windows  are  inviolable, 
whereas  soft  windows  may  be  violated  at  a  cost. 

For  both  hard  time  windows  and  soft  time  windows,  we 
developed  some  new  heuristics,  SCCC  and  SCAO,  which  are 
modifications  of  Stewart's  unconstrained  TSP  heuristics 
[Ref.  6]  CCCO  and  CCAO.  Also  for  the  hard  time  window  only, 
we  developed  the  SLACK  heuristic.  We  also  developed  an 
exact  algorithm  for  both  hard  and  soft  window  using  state 
space  relaxation  dynamic  programming  and  branch  and  bound  as 
proposed  by  Chr istof ides  et  al.  [Ref.  .5]. 

The  procedures  were  shown  to  be  effective  on  some  moder- 
ately small  sized  problems.  A  Nearest  Neighbor  heuristic 
was  also  developed,  cut  it  was  often  unable  to  solve  the 
problem  with  hard  time  windows,  and  it  found  very  low 
guaiity  solutions  with  soft  time  windows.  Tnis  experience 
is  consistent  with  tne  findings  of  ethers  [Ref.  7]  who 
determined  that  the  Nearest  Neighbor  heuristic  does  not 
perform  well  en  the  unconstrained  TSP. 

The  SCCO  and  SCAO  are  generally  effective  en  most  of  the 
small  sized  problems  we  tested,  except  for  tne  problems  in 
which  less  than  half  the  nodes  have  time  windows-  Further 
research  is  needed  in  order  to  satisfactorily  solve  these 
problems.  Another  problem  difficulty  that  may  reguire  more 
research  is  dealing  with  wider  time  windows. 

The  SLACK  heuristic  which  is  used  only  with  hard  time 
windows   is   slightly   slower   than   the   otner   heuristics. 
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Particularly,   in  the  lower  time  window  percentage  problems, 
the  accuracy  becomes  lower. 

The  exact  algorithm  succeeded  in  solving  1 4  of  the  16 
test  problems  to  optimality,  but  it  was  too  slow  to  use  in 
most  of  the  lower  time  window  percentage  problems.  This 
algorithm's  performance  also  depends  upon  the  quality  of  the 
upper  bound  which  is  obtained  from  the  heuristic. 
Additional  research  is  needed  to  reduce  computation  time, 
but  •  a  working  program  for  at  least  some  problems  has 
resulted  from  this  effort. 
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Depot    cc-ordinates   :     (12,10) 

problem   source  :    Sedgewick   [fief.    19] 
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APPEBDIX    B 
TEST   PROBLEM    [2] 
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Depot   co-ordinates   :     (326,181) 

problem   source    :     Eilon  et   al-    [Ref.    21] 


80 


APPEBDIX    C 
TEST   PEOBLEH    [3] 
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Depot   co-ordinates   :     (145,215) 

problem   source   :    Eilon  et  al.    [fief.    21], 
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APPJjJII   D 

TEST   PROBLEM    £5] 
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Depot  co-ordinates   :     (30,40) 

problem    source    :    Ellon  et  al.    fRef    .21]. 
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APPE8DIX    E 
TEST   PBOBLEH   £6] 
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Depot    co-ordinates  ;      (40,40) 
problem  source    :    Eilon   et  al.    [ fief   .21]. 
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APPEHDIX    F 
TEST    PfiOBLEM    FOE   THE    SCCO 
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Depot    co-ordinates   :      (12,10) 

problem   source 

node    locations   :     Sedgevick    [Ref.     15] 
time      windows      :    see    Chapter    V- 
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i££EBDIX  G 
TEST  PROBLEM  [ 1-1 ] 
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17 

2    11 

1 

46 

53 

I   12 

16 

14 

2 

9 

3     6 

8 

27 

36 

I   13 

15 

2 

51 

59 

4     4 

3 

37 

45 

I    14 

13 

16 

5 

13 

5     5 

15 

18 

28 

1   15 

2 

12 

22 

30 

6     8 

11 

14 

23 

I   16 

12 

10 

- 

- 

7     6 

4 

35 

43 

8     7 

4 

42 

49 

9     9 

7 

58 

68     | 

0    14 

5 

53 

64 

Depot  co-ordinates  :  (12,10) 
CL  =  2.  0,  CU  =5.0 
problem  source 

node  locations  :  Sedgewick  [Ref. 

tine   windows   :  see  Chapter  V. 


19] 
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APPEBDIX    H 
TEST    PROBLEM   [ 1-2  ] 


|  node 

X 

y  time 

window   , 

node 

X 

y 

time 

window   J 

l(i> 

u(i) 

Ki) 

a(i)   1 

I    1 

3 

9 

25 

32 

11 

10 

13 

10 

17    I 

I    2 

11 

1 

46 

53 

I   12 

16 

14 

2 

9    I 

I    3 

6 

8 

27 

36 

I   13 

15 

2 

51 

59    | 

I    4 

4 

3 

- 

- 

I   14 

13 

16 

5 

13    | 

1    5 

5 

15 

18 

28     | 

15 

2 

12 

22 

30         j 

1    6 

8 

11 

14 

23 

16 

12 

10 

- 

I 

I    7 

6 

4 

35 

43 

1    3 

7 

4 

42 

49 

1    9 

9 

7 

58 

68 

1   10 

14 

5 

53 

64 

Depot  co-ordinates   :     (12,10) 

CL    =    2.  0,    CU    =5.0 

problem  source 

node  locations  :  Sedgewicx  [Ref 
time   windows   :  see  Chapter  V. 


19] 
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APPEJDIX    I 
TEST    PEOBLEH   [ 1-3 ] 


node 

X 

y  time 

window 

node 

X 

y 

time 

window   | 

Mi) 

u(i)    | 

Mi) 

Mi)   J 

1 

3 

9 

— 

- 

11 

10 

13 

10 

17    J 

2 

11 

1 

46 

53 

12 

16 

14 

2 

9    I 

3 

6 

8 

27 

36 

I   13 

15 

2 

51 

59    I 

4 

4 

3 

37 

45 

I   14 

13 

16 

5 

13    l 

5 

5 

15 

18 

28     | 

15 

2 

\2 

22 

20         I 

6 

8 

11 

14 

23     | 

16 

12 

10 

- 

j 

7 

b 

4 

35 

43 

8 

7 

4 

42 

49     , 

9 

9 

7 

- 

- 

10 

14 

5 

- 

- 

Depot  co-ordinates   :     (12,10) 
CL    =    2-  0,    CU    =5.0 
problem   source 

node    locations    :     Sedgewick    [Ref. 
tine      windows      :    see   Chapter    V. 


19] 


87 


A£P EHDIX    J 
TEST    PEOBXEH   [ 1-4  ] 


I    node 

X 

y  time 

window 

node 

X 

y 

time 

window   i 

3 

KiJ 

u(i)    J 
32 

I   11 

10 

13 

Mi) 

10 

a(i)   I 

|    1 

9 

25 

17    i 

I    2 

11 

1 

46 

53 

I   12 

16 

14 

- 

j 

J    3 

0 

8 

- 

- 

I   13 

15 

2 

51 

59    J 

I    4 

4 

3 

37 

45 

I   14 

13 

16 

5 

13    | 

J    5 

5 

15 

- 

- 

15 

2 

12 

- 

-    | 

I    6 

8 

11 

14 

23 

I   16 

12 

10 

- 

i 

|    7 

6 

4 

- 

- 

I         8 

7 

4 

- 

- 

I    9 

9 

7 

58 

68 

I   10 

14 

5 

- 

- 

Depot   Co-ordinates   :     (12,10) 
CL    =    2-  0,    CU    =5.0 
Problem    Source 

node   locations    :    Sedgewick   (Ref. 

time    windows         :   see    Chapter    V. 


19] 


88 


APPEBDII    K 
TEST    PROBLEM    [2-1  ] 


J  node 

X 

y 

tide  window 

node 

X 

y   time  window 

27  2 

Mi) 
125 

Mi)    I 
135    | 

12 

267 

242 

Mi) 
170 

u(i) 

|    1 

295 

179 

I    2 

301 

25  8 

110 

118 

I   13 

259 

265 

193 

202 

I    3 

309 

260 

102 

1  10    j 

I   14 

315 

2  33 

57 

67 

I    4 

217 

274 

242 

250 

I   15 

3  29 

252 

81 

89 

|    5 

218 

27  8 

23  9 

246 

I   16 

318 

252 

90 

98 

I    6 

282 

26  7 

141 

149 

I   17 

3  29 

224 

40 

49 

I    7 

242 

249 

279 

2  86 

I   18 

2  67 

213 

382 

393 

I    8 

230 

26  2 

261 

271 

I   19 

275 

192 

404 

413 

I    9 

249 

26  8 

206 

215 

20 

3  03 

201 

432 

442 

I   10 

256 

26  7 

200 

208 

I   21 

208 

217 

3  23 

332 

I   11 

265 

257 

183 

193 

|   22 

3  26 

181 

- 

- 

Depot   co-ordinates   :     (326,181) 
CL    =    2.  0,    CU    =5.0 
problem   source 

node    locations    :     Eilon   et    al-    [  Ref .    21], 
time      windows      :     see   Chapter   V. 
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APPEBDIX    L 
TEST    PROBLEM    [2-2] 


J    node 

X 

y 

time    window 

node 

X 

y      time    window      J 

KiJ 

Mi)          I 

Ki) 

a  (i )       i 

1         1 

295 

112 

125 

135         | 

12 

267 

242 

170 

179         j 

I         2 

301 

25  8 

110 

118 

I       13 

259 

265 

- 

j 

|         3 

309 

26  0 

102 

1  10 

I       14 

315 

233 

57 

67          | 

i             4 

217 

27  4 

242 

250 

I       15 

3  29 

2  52 

81 

89          | 

I         5 

218 

27  8 

23  9 

246 

I        16 

318 

2  52 

90 

98          l 

I         6 

282 

267 

141 

149 

I       17 

329 

224 

40 

49          | 

I         7 

242 

249 

279 

286 

I        18 

2  67 

213 

382 

393          | 

I         8 

230 

26  2 

261 

271 

I       19 

275 

192 

404 

413          ] 

I         9 

249 

26  8 

206 

215         | 

20 

3G3 

201 

432 

442          j 

|        10 

25b 

26  7 

20  0 

208 

I       21 

208 

217 

- 

| 

i        11 

265 

25  7 

183 

193 

|       22 

326 

181 

- 

J 

Depot    co-ordinates   :     (326,181) 

CL    =    2-  0,    CU    =5.0 

problem   source 

node    locations    :     Eilon   et    al.    [Ref.    21] 
tiioe      windows      :     see   Chapter   V    . 
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APPEBDIX    H 
TEST    PBOBIEH   [2-3  ] 


|  node 

X 

y 

time  window 

node 

X 

y   time  window   j 

Ki) 

u<i) 

Ki) 

u(ij   J 

|    1 

295 

272 

125 

135    | 

12 

267 

242 

170 

179    1 

I    2 

301 

25  3 

- 

- 

I   13 

259 

265 

- 

{ 

I    3 

309 

26  0 

- 

- 

I   1<* 

315 

2  33 

57 

67    | 

I    4 

217 

274 

242 

250 

I   15 

3  29 

252 

81 

89    | 

I    5 

218 

278 

239 

246 

I   16 

318 

252 

90 

95    1 

I    6 

282 

26  7 

141 

149 

I   17 

3  29 

224 

40 

49    j 

I    7 

242 

24  9 

- 

- 

I   18 

267 

213 

382 

393    1 

I    3 

230 

262 

- 

I 

19 

275 

192 

- 

j 

I    9 

249 

26  8 

206 

215    | 

20 

3  03 

201 

432 

442    1 

I   10 

256 

26  7 

200 

208 

I   21 

208 

217 

3  23 

332    J 

I   11 

265 

257 

183 

193 

22 

326 

181 

- 

-     j 

Depot   co-ordinates   :     (326,181) 
CL    =    2-  0,    CU    =   5.0 
problem   source 

node    locations    :    Eilon   et    al.    [Ref.    21]. 

tine      windows      :     see  Chapter   V. 
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APP E8DIX    N 
TEST    PBOBLEM   [2-4  ] 


|  node 

X 

y 

time  window 

node 

X 

y   time  window   j 

Ki) 

u(i) 

Mi) 

u(i)   j 

j    1 

295 

272 

- 

- 

12 

267 

2  42 

170 

179    j 

I    2 

301 

25  8 

110 

1  18 

I   13 

259 

265 

- 

j 

|    3 

309 

26  0 

- 

- 

I   14 

315 

233 

57 

67    J 

<    4 

217 

274 

242 

2  50 

I   15 

329 

252 

- 

| 

i          5 

218 

27  8 

239 

246 

16 

318 

252 

90 

98    1 

I     o 

282 

26  7 

141 

149 

I   17 

329 

224 

- 

j 

i    7 

242 

24  9 

279 

286 

I   18 

267 

213 

- 

j 

I    8 

23  0 

26  2 

- 

- 

I   19 

275 

192 

404 

413    J 

I    9 

249 

268 

206 

215 

20 

3  03 

201 

432 

442    | 

|   10 

256 

26  7 

- 

- 

I   21 

208 

217 

- 

j 

I   H 

265 

257 

- 

- 

|   22 

3  26 

181 

- 

| 

Depot  cc-orainates 
CL  =  2.  Q,  CU  =5.0 
problem   source 

node   locations    : 
time      windows      : 


(326,  181) 


£ilon  et  al.  [  Ref .  21  ]. 
see  Chapter  V. 
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APP E1DIX    0 
TEST    PROBLEM   [3-1  ] 


node 

X 

y 

time  window 

node 

X 

y   time  window   j 

Ki) 

Mi)    I 

1U) 

a(i)   i 

1 

151 

264 

196 

204    | 

12 

156 

217 

105 

118    j 

2 

159 

261 

185 

193 

I   13 

129 

214 

259 

271    | 

3 

130 

254 

217 

2  25    | 

14 

146 

2  08 

2 

10    J 

4 

128 

252 

222 

234 

I   15 

164 

208 

92 

105    I 

5 

163 

247 

174 

185    | 

16 

141 

206 

10 

19    | 

6 

146 

24  6 

142 

154 

17 

147 

193 

54 

68    | 

7 

161 

242 

166 

173 

|   18 

164 

193 

79 

89    J 

8 

142 

23  9 

131 

142 

I   19 

129 

139 

30 

38    | 

9 

163 

236 

159 

165    | 

20 

155 

185 

67 

75    J 

10 

148 

23  2 

123 

131 

I   21 

139 

182 

^0 

53    | 

11 

128 

23  1 

242 

2  53 

22 

145 

215 

- 

j 

Depot    co-ordinates   :     (145,2  15) 
CL    =    2.  0    ,    CO    =    5-0 
problem   source 

ncde    locations    :    Eilon   et    al.    [Eef 
time      windows      :    see  Chapter   V. 


21]. 
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APPE1DIX    P 
TEST    PROBLEH   [3-2  ] 


|  node 

X 

y 

time  window   , 

node 

X 

y   time  window   | 

l(i) 

Mi)    I 

Mi) 

a  (i )   I 

|    1 

151 

264 

196 

2  04 

12 

156 

217 

105 

118    | 

I    2 

159 

26  1 

185 

193 

I-   13 

129 

214 

- 

| 

|    3 

130 

254 

217 

225 

14 

146 

2  08 

2 

10    | 

I    4 

128 

25  2 

222 

234 

I   15 

164 

208 

92 

105    | 

|    5 

163 

247 

174 

185    | 

16 

141 

206 

10 

19    j 

I    6 

146 

24  b 

142 

154    | 

17 

147 

193 

54 

68    J 

I    7 

161 

242 

166 

173 

|   18 

164 

193 

79 

89    1 

I    8 

142 

23  9 

131 

142 

I   19 

129 

189 

30 

38    | 

I    9 

163 

236 

159 

165    | 

20 

155 

185 

67 

75    J 

I   10 

148 

23  2 

123 

131 

I   21 

139 

182 

- 

| 

i   11 

128 

23  1 

242 

253 

22 

145 

215 

- 

j 

Depot    co-ordinates   :     (145,2  15) 
CL    =    2.0    ,    CU    =    5-0 
problem    source 

node    locations    :     Eilon   et    al   [Ret.    21]. 
titte      windows      :     see   Chapter   V. 


94 


TEST    PROBLEM   £3-3  ] 


node 


1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 


X 

y 

time  window 

node 

X 

y   time  window   j 

i(i) 

u(i)    I 

Mi) 

u(i)   I 

151 

264 

196 

204    | 

12 

156 

217 

105 

118    | 

159 

26  1 

- 

- 

I   13 

129 

214 

- 

j 

130 

254 

- 

| 

14 

146 

2  08 

2 

10    | 

123 

25  2 

222 

234 

I   15 

164 

208 

92 

105    | 

163 

247 

174 

185    | 

16 

141 

2  06 

10 

19    j 

146 

24  6 

142 

154 

17 

147 

193 

54 

68    | 

161 

242 

- 

- 

18 

164 

193 

79 

89    j 

142 

23  9 

- 

- 

19 

129 

189 

- 

j 

163 

236 

159 

165    | 

20 

155 

185 

67 

75    | 

148 

23  2 

123 

131 

I   21 

139 

182 

40 

53    | 

128 

23  1 

242 

253 

22 

145 

215 

- 

| 

Depot    co-ordinates   :     (145,2  15) 
CL    =    2. 0    ,    CU    =    5.0 
problem   source 

code    locations    ;    Eilon   et    al.    [Ref.    21], 
time      windows      :     see  Chapter  V. 
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APPE8DIX    H 
TEST    PEOBLBH   £3-4  ] 


node 

X 

y 

time  window 

node 

X 

y   time  window   J 

Mi) 

u(i) 

Mi) 

u(i)   1 

1 

151 

264 

- 

- 

12 

156 

217 

105 

113    | 

2 

159 

26  1 

.  185 

193 

I   13 

129 

214 

- 

j 

3 

130 

254 

- 

- 

14 

146 

2  08 

2 

10    | 

4 

128 

25  2 

222 

234 

1   15 

164 

208 

- 

j 

5 

163 

24  7 

174 

185 

16 

141 

206 

10 

19    | 

6 

146 

24  6 

142 

154 

17 

147 

193 

- 

j 

7 

161 

242 

166 

173 

18 

164 

193 

- 

| 

8 

142 

23  9 

- 

- 

19 

129 

189 

30 

38    | 

9 

163 

236 

159 

165 

20 

155 

135 

67 

75    | 

10 

148 

23  2 

- 

- 

I   21 

139 

182 

- 

| 

11 

128 

23  1 

- 

- 

|  22 

145 

215 

- 

j 

Depot   co-ordinates   :     (145,215) 

CL    =    2. 0    ,    CU    =    5.0 

problem   source 

node    locations    ;    Eilon    et    al.    [Ret.    21], 
time      windows      :     see  Chapter   V. 
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APPEBDIX   S 
TEST    PROBLEM   £4-1  ] 


|  node 

X 

y 

time  window   | 

node 

X 

y   time  window   l 

Mi) 

u(i) 

Mi) 

Mi)   J 

I    1 

151 

264 

171 

179 

12 

156 

217 

72 

79    | 

I    2 

159 

261 

162 

170 

I   13 

129 

214 

237 

245    j 

1    3 

130 

254 

196 

203 

14 

146 

2  08 

5 

9    l 

1    4 

128 

252 

198 

206 

15 

164 

208 

61 

67    j 

|    5 

163 

24  7 

128 

136    | 

16 

141 

2  06 

10 

14    j 

1    6 

146 

24  6 

106 

1  13 

17 

147 

193 

22 

28    j 

|    7 

161 

242 

122 

130 

|   18 

164 

193 

48 

53    | 

1    3 

142 

23  9 

97 

105 

I   19 

129 

189 

261 

2b9    | 

1    9 

163 

236 

138 

146    j 

20 

155 

185 

35 

40    J 

1   10 

148 

232 

89 

96 

i   21 

139 

182 

273 

280    | 

I   11 

128 

23  1 

220 

2  27 

|   22 

145 

215 

- 

| 

Depot    co-ordinates   :     (145,215) 

CL    =    2. 0r      CU    =    5-0 

problem    source 

node    locations    ;    Zilon    et    al   [  fief .    21], 
time    windows        :     see  Chapter   V. 
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APP JJDIX    T 
TEST    PfiOBlEH   [4-2] 


1    node 

X. 

y 

time    window 

node 

X 

y      time    window 

l(i) 

u(i)          I 

Mi) 

u(i) 

J         1 

151 

264 

171 

179 

12 

1  56 

217 

72 

79 

1         2 

159 

26  1 

162 

170 

13 

129 

214 

- 

- 

1         3 

130 

254 

196 

203         | 

14 

146 

2  08 

5 

9 

1           4 

128 

252 

196 

2  06 

15 

164 

208 

61 

57 

|         5 

163 

24  7 

128 

136 

16 

141 

2  06 

10 

14 

1         t> 

146 

24  6 

10b 

113          | 

17 

147 

193 

22 

28 

1          7 

161 

242 

122 

130 

I        18 

164 

193 

48 

53 

1         8 

142 

23  9 

97 

105 

I        19 

129 

189 

261 

269 

J          9 

163 

236 

138 

146          j 

20 

155 

185 

35 

40 

1       10 

148 

232 

69 

96 

I       21 

139 

182 

- 

- 

I        11 

128 

231 

220 

2  27 

22 

145 

215 

- 

- 

Depot   co-ordinates   :     (145,2  15) 

CL    =    2.0,       CU    =    5.0 

problem    source 

node    locations    :    Eilon   et    al.    [3ef.    21], 
time      windows      :     see   Chapter    V. 
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APPEMDIX  0 
TEST  PBOBLEH  [4-3  ] 


node 


1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 


X 

y 

time    window 

node 

X 

y      time    window      J 

Mi) 

u(i) 

Ki) 

u(i)       I 

151 

26a 

171 

179 

I       12 

156 

217 

72 

79          | 

159 

261 

- 

- 

I       13 

129 

214 

- 

| 

130 

254 

- 

- 

I       1" 

146 

2  08 

5 

9          I 

128 

252 

198 

206 

1       15 

164 

208 

61 

67          | 

163 

247 

128 

136 

16 

141 

206 

10 

14          | 

146 

246 

106 

1  13          j 

I       17 

147 

193 

22 

28          | 

161 

242 

- 

- 

I       18 

164 

193 

48 

53          | 

142 

23  9 

- 

- 

I       1^ 

129  • 

189 

- 

j 

163 

236 

138 

146          | 

20 

155 

185 

35 

40          | 

148 

23  2 

89 

96 

I       21 

139 

182 

273 

280          | 

128 

231 

220 

2  27 

I       22 

145 

215 

- 

j 

Depot  cc-ordinates  :  (145,2  15) 

CL  =  2.0,       CO  =  5.0 

problem    source 

node    locations    :    Eilon   et   al.    [Ref.    21] 
time    windows        :     see  Chapter   V. 
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1ZEJJDIX    ¥ 
TEST    PROBLEM   [4-4 ] 


j  node 

X 

y 

time  window   | 

node 

X 

y   time  window 

Mi) 

tt'(i)    I 

Mi) 

u(i) 

1    1 

151 

264 

- 

- 

12 

156 

217 

72 

79 

1    2 

159 

261 

1c2 

170 

13 

129 

214 

- 

- 

1    3 

130 

254 

- 

- 

14 

146 

2  08 

5 

9 

1    4 

128 

25  2 

198 

2  06 

I   15 

164 

208 

- 

- 

1    5 

163 

247 

128 

136 

16 

141 

2  06 

10 

14 

1    6 

146 

24  6 

106 

113 

17 

147 

193 

- 

- 

1    1 

161 

242 

122 

130 

I   18 

164 

193 

- 

- 

1    3 

142 

23  9 

- 

- 

I   19 

129 

189 

261 

269 

1    9 

163 

236 

138 

146    | 

20 

155 

185 

35 

40 

1   10 

148 

23  2 

- 

- 

I   21 

139 

182 

- 

- 

1   11 

128 

231 

- 

- 

22 

145 

215 

- 

- 

Depot    co-ordinates   :     (1<*5r215) 
CL    =    2. 0,       CU    =    5.0 
problem    source 

node    locations    :    Eilon    et    al.    [fief.    21]. 

time      windows      :     see   Chapter   V. 
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